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ABSTRACT 


Two-phase turbines open the possibility of new types of power cycles 
operating with extremely wet mixtures of steam and water, organic fluids, or 
immiscible liquids and gases. Possible applications are geothermal power, 
waste-heat recovery, refrigerant expansion, solar conversion, transportation 
turbine engines, and engine bottoming cycles. 

A theoretical model for two-phase impulse turbines was developed. 

Apparatus was constructed for testing one- and two-stage turbines (using speed 
decrease from stage to stage). Turbines were tested with water-and-ni trogen 
mixtures and Refrigerant 22. Nozzle efficiencies were 0.78 (measured) and 
0.72 (theoretical) for water-and-ni t rogen mixtures at a water/nitrogen mixture 
ratio of 68, by mass; and 0.89 (measured) and 0.84 (theoretical) for Refrig- 
erant 22 expanding from 0.02 quality to 0.28 quality. Blade efficiencies 
(shaft power before windage and bearing loss divided by nozzle jet power) were 
0.63 (measured) and 0.71 (theoretical) for water-and-ni trogen mixtures and 0.62 
(measured) and 0.63 (theoretical) for Refrigerant 22 with a single-stage 
turbine, and 0.70 (measured) and 0.85 (theoretical) for water-and-ni trogen 
mixtures with a two-stage turbine. 



PRECEDING PAGE DLANX NOT ptl-MED 


ACKNOWLEDGMENT 

The two-phase turbine research was supported from May 1977 to September 
1980 by the De;*artment of Energy (DOE), Division of Geothermal Energy, Hydro- 
thermal Technology Branch, Clifton B. McFarland, Chief, uncjr Interagency 
Agreement EG-78-A-30-1026 (NASA Task Order RD 152, Amendment 130), “Analysis 
Support in Geothermal Cycles and Expanders." The turbine research continued 
with support from the DOE Division of Advanced Energy Projects, Dr. Ryszard 
Gajewski, Director, under Interagency Agreement DE-AI01-80ER1061A (NASA Task 
Order RD 152, Amendment 266), "Two-Phase Turbines for Efficient Waste Recovery." 

This report is the final repoit on the work for the Division of Geothermal 
Energy and the midterm report on the work for the Division of Advanced Energy 
Pro jec t s. 

The turbine research was conducted in the Acoustics, Fluid Dynamics, and 
Heat Transfer Group (Paul F. Massier, Supervisor) of the Energy and Materials 
Research Section (Dr. Robert F. Landel, Manager). Support in organizing the 
project and coordinating with NASA and DOE was provided by G.W. Meisenholder, 
Manager of Transportation Programs at JPL. Special acknowledgment is due P.R. 
Miller, Manager of Conservation and Fossil Energy Systems at NASA, for assis- 
tance in securing support for the research and handling contractual arrange- 
ments. 


V 


PRECEDING PAGE CLANK NOT PRI MED 


CONTENTS 

I. INTRODUCTION 1 

II. BACKGROUND 2 

III. APPLICATIONS 4 

A. Saturated Water Turbines 4 

B. Organic-Fluid Turbines 4 

C. Wet-'to-Dry Cycle 12 

D. Two-Component Cycle 15 

IV. NOZZLES 19 

A. Construction 19 

B. Nozzle Program 22 

C. Performance Definitions 24 

D. LLL Steam-and-Water Data 26 

E. Wa ter-and-Nit rogen Data 27 

F. Refrigerant-22 Data 32 

G. LLL Turbine Nozzle 32 

H. Wet-To-Dry Nozzle 37 

V. ROTORS 39 

A. Erosion 39 

B. Liquid Prth 40 

C. Rotor Program 4 7 

D. Rotor Test Apparatus 31 

E. Measurement Procedure 34 

F. Water-and-Nit rogen Turbine Data 54 

G. Refrigerant-22 Turbine Data 56 

H. LLL Steam-and-Water Turbine Data 60 

I. Staging Method 61 

J. First-Stage Rotor Data b5 

vi i 


K* Two*-Stage Turbine Data * * • • 68 

L« Divergence and Stagnation Losses * • 79 

M* Separator Turbine Concept ••••«• ***«•*« 85 

N. Turbine Performance Examples. * ... * 89 

VI. CONCLUSION 96 

A. Analysis Capability 96 

B. Application Prospects 96 

APPENDIXES 

A. ROTOR MODEL 99 

B. SEPARATOR TURBINE MODEL I33 

REPERENCES 139 

Figures 

1. Steamrand^water cycle 5 

2. Organic-fluid two-phase turbine ... 6 

3. Organic-fluid two-phase cycle . 8 

4. Comparison of Rankine and t^/o-phase cycles 10 

5. Comparison of supercritical and two-phase cycles . 11 

6. Steam topping cycle with organic-fluid bottoming cycle. .... 13 

7. Wet-to-dry cycle ..... «..•••.••«•• 14 

8. Two-phase, two-component engine . 16 

9. Two-phase, two-component cycle .... ..... 17 

10. Two-phase nozzle . 20 

11* Experimental two-phase nozzle *..••• 21 

12. Water-and-nit rogen jet at 80 m/s 23 

13. Comparison of theoretical and experimental drop diameters 

and velocity coefficients for the LLL steam-and-water nozzle • 28 

14. Comparison of theoretical and experimental exit velocities . . 29 

15* Wet-to-dry expansion of toluene from 289^C to 70^C .... 38 

16. Effect of drop diameter on liquid impingement 41 

viii 


17. Rotor 1 blade ahape 42 

18. Liquid exit flow at high liquid/gaa mixture ratio . . 43 

19. Axial view of exit flow atR«45 .......... 45 

20. Exit flow atR*>ll.... 46 

21. Exit flow atR*>4 48 

22. Two-phase turbine experiment .... 52 

23. Single-stage turbine 53 

24. Comparison of theoretical and experimental torques of 

water-and-nitrogen turbine . 55 

25. Comparison of theoretical and experimental efficiencies 

of water-and-nitrogen turbine 57 

26. Comparison of theoretical and experimental torques of 

LLL steam-and-water turbine ......... .. 62 

27. Comparison of theoretical and experimental efficiencies 

of LLL steam-and-water turbine 63 

28. Staging method 66 

29. Rotor 2 blade shape ... ............ 67 

30. Exit flow from Rotor 2 69 

31. Axial view of exit flow from Rotor 2 70 

32. Comparison of exit flow from flat and curved blades ...... 71 

33. Comparison of theoretical and experimental blade torques 

of Rotor 2 72 

34. Two-stage turbine .... ..... 73 

35. Effect of nozzle position on second-stage zero-speed 

torque at 2200-rpm first-stage speed ............. 75 

36. Comparison of theoretical and experimental second-stage 

torques of two-stage turbine at 2200-rpm first-stage speed . . 76 

37. Comparison of theoretical and experimental efficiencies of 

two-stage turbine at 2200-rpm first-stage speed ........ 77 

38. Exit flow from a single blade 80 

39. Exit jet and delayed drainage flow from Rotor 2. 82 

40. Exit flow from bladed arm ................... 83 

41. Rotor 2 exit scoop for flow measurement 84 

42. Separator turbine 86 


A3* Conparison of efficiencies of impulse turbine end 

sepereCor turbine .... 88 

AA* Theoretical efficiencies of S-MW steam-and>water and 

Refrigerant-113 turbines .. ....... 95 

A-1. Two-phase turbine geometry . lOO 

A-2. Definitions of velocity vectors .......... 102 

A-3. Impingement geometry ...... ....... 104 

A-A. Impact behavior ...... io6 

A-5. Impact and friction on a blade section io8 

A-6. Blade lengths ....... 113 

B-1. Separator turbine 13A 


Tables 


1. Comparison of theoretical and experimental nozzle 

performance using water and nitrogen ....... .30 

2. Comparison of theoretical and experimental nozzle 

performance using Refrigerant 22 33 

3. Theoretical performance of LLL turbine nozzle ......... 35 

A. Comparison of theoretical and experimental turbine 

performance using water and nitrogen 58 

5. Comparison of theoretical and experimental turbine 

performance using Refrigerant 22 .... 59 

6. Comparison of theoretical and experimental performance 

of LLL steam-and-water turbine 6A 

7. Comparison of theoretical and experimental two-stage 

turbine performance . ..... 78 

8. Theoretical performance of 5-MW steam-and-water turbine .... 90 

9. Theoretical performance of 5-MW Refrigerant-113 turbine .... 93 

A-1* Rotor model nomenclature ... .....125 


X 


I. INTRODUCTION 


Power conversion using two-phase liquid-gas mixtures has been investigated 
for use in pumping (Ref* 1), liquid-metal magnetohydr.>dynamics (Ref. 2), and 
turbine engines (Ref. 3). ^is report describes current research on two-phase 
turbines at JPL. 

Two-phase turbines provide a way to operate in the wet region of water-and- 
steam mixtures or of organic fluids. In addition) two-phase, two-component 
mixtures such as steam and oil or steam and molten salts can, in principle, be 
used as working fluids in unusual cycles. Possible applications of two-phase 
turbines are in geothermal power, waste-heat recovery, refrigerant expansion, 
solar conversion, transportation turbine engines, and engine bottoming cycles. 
Potential advantages over conventional vapor cycles in these applications are 
more effective use of the energy in such hot liquids as geothermal fluids, 
better matching to sensible-heat sources, lower shaft speeds, and the 
flexibility of operating without the restriction of dry-vapor expansion. 

This report will present theory and test results for two-phase turbines 
and give examples of applications. 
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II. BACKGROUND 


Earlier, I reported two-phese flow research on jet punpa (Ref. 1) and 

liquid~aietal nagnetohydrodynamics (Ref. 2). Those two projects were based on 

the ability of a two-phase flow to provide pumping without using moving 

mechanical parts. By expanding a two-phase mixture in a noxsle and then 

removing the gas phase, a large rise in the dynamic pressure of the fluid can 

be achieved. The separation of gas from liquid increases the density terms 

2 

in the dynamic pressure ijpV , enabling the separated liquid to flow 
through a diffuser to reach a pressure much higher than the inlet pressure of 
the two-phase nozzle. 

Instead of a pressure rise, power can be taken out of the liquid. This 
effect is the basis of liquid-metal magnetohydrodynamics, or liquid-metal 
MHD. A liquid metal is circulated in a closed loop, and power is extracted 
from the liquid metal with a magnetic field. Liquid-metal MHD has been studied 
at the Jet Propulsion Laboratory (JPL), the Argonne Nationsl Laboratory, and 
elsewhere. Because of the efficiency limitations of two-phase nozzles, 
separators, diffusers, and generators, liquid-iiietal MHD cycles have been 
limited to about half the efficiency of turbine cycles, and applications 
appear limited to extreme temperature conditions where turbines cannot be used. 

During the JPL liquid-metal MHD project in 1973, I suggested, not 
seriously, that the MHD generator could be replaced by a two-phase liquid- 
metal turbine to raise the efficiency of the liquid-metal MHD cycle. Lance 
Hays, an engineer on the project, recognized in this suggestion the germ of a 
serious possibility and pointed out that a two-phase turbine, coupled with a 
pump to return the liquid to the nozzle, could use two-phase mixtures other 
than liquid metals, such mixtures as steam and hydraulic oil, to provide a new 
type of turbine engine with low shaft speed and other possible «dvantages. 

At the same time, A.L. Austin and his geothermal group at Lawrence 
Livermore Laboratory (LLL) proposed a two-phase turbine for geothermal power 
(Ref. A). The two-phase fluid from geothermal wells (both the water and the 
steam) would be passed through the turbine. The LLL group pursued this idea 
experimentally in research on two-phase impulse turbines (Ref. 5). They 


hoped thee the veter drop eites would be emell enough for the eteem end weter 
to beheve ee a homogeneous fluid end give these turbines efficiencies 
compereble to *«hose of vepor turbines* 

Besed on celculetions of drop sites end drop trejectorieSi I concluded 
thet the liquid drops in tv -phese turbines would be too lerge to follow the 
ges phese, end thet the liquid would impinge on the bledes end form e thin 
liquid film* The friction dreg of this film would be so high thet the 
efficiency of single*stege impulse turbines would be limited to little more 
then 30 percent* 

To Achieve higher efficiency , Keys end I thought it would be necessary to 
use a '^seperetor turbine^* in which the liquid would be sepereted from the 
vepor in e rotary separator and then passed through a liquid turbine* This 
concept was subsequently pursued by Hays' company . Biphase Energy Systems, in 
seeking to develop comnercial applications (Refs* 6 end 7}* 

In early turbine experiments at JPL during 1978, I studied the flow of 
water-and*»nitrogen mixtures in rotary separators and concluded thet the impact and 
other losses in rotary separators would limit turbine efficiency as severely 
as the liquid friction in bladed turbines* At the same time, LLL furnished 
JPL a surplus steam turbine for bladed turbine tests* The second row of 
blades was cut away to make a single-stage rotor* In water-and-nitrogen tests, 
this rotor gave as high an efficiency as seemed attainable with a separator 
turbine* In addition, I found it possible to further improve the efficiency 
by using two stages, with high speed in the first stage and a lower speed in 
the second stage* This two-phase turbine is the type discussed in this report* 

It is also possible to achieve efficient two-phase expansion using 
positive-displacement expanders* Helical'^screw (Lysholm) expanders have been 
tested by LLL (Ref* 8) using a two-phase steam-and-water flow, and by Hydrothermal 
Power Company in cooperation with JPL (Ref* />* Efficiencies as high as 60 
percent have been achieved in geothermal field tests at 300-kW output power 
(Ref* 9)* 
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III. APPLICATIOHS 


A. Saturated Water TUrbinee 

One of the siapleet applicatione of two-phaae turbinea ia the open^ycle 
expanaioo of aaturated water or lo%r-quality ateam (^'quality" ia the ateaai 
fraction by maaa). An example ia the expanaioo of fluid from a geothermal 
well. 


Figure 1 ahowa ateam-and~water expanaioo on a temperature'-entropy (T~’S) 
diagram. The expanaioo atarta on or near the aaturated liquid line and 
proceeda doiraward to the ateam condenaer preaaure. Thia ia the proceaa 
propoaed aa the "total^flow" concept by LLL (Refa. 4 and S). 

An alternative approach for geothermal uae ia to replace with two-pbaae 
turbinea the throttling atepa uaed to provide ateam in the conventional 
geothermal cycle. The output from the two-phaae turbinea adda to the output 
of the ateam turbinea. Haya and CerinI (Refa. 6 and 7) have developed auch 
geothermal expandera, combining flaah vaporization, aeparation. and power 
output in a aingle unit. 

B. Organic-Fluid Turbinea 

One of the moat promiaing applicationa of two-phaae turbinea ia in 
cloaed-cycle enginea uaing organic working fluida. The advantagea over 
Rankine vapor cyclea are better matching to the cooling curve of aenaible-beat 
aourcea and elimination of the boiler. Poaaible applicationa arc geothermal 
cloaed-cycle (binary) power planta, engine cxhauat heat recovery, induatrial 
waate-faeat recovery, and bottoming cyclea for ateam and gaa turbine planta. 

Figure 2 ahowa the equipment arrangement in an organic-fluid two-pbaae 
turbine engine. The fluid to be cooled (the *'aource fluid*'), auch aa 
geothermal hot water or engine exhauat. flowa through the heat exchanger from 
A to B and ia cooled by counterflow heat tranafer to the organic %rarking 
fluid. The working fluid ia hea'.ed to saturation temperature. The aaturated 
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Fig. 2. Organic-fluiu two-phase turbine 
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liquid, with possibly a small amount of vapor, flows to a two~phase nocala. 

The liquid expands, partially vaporises, and accelerates in the nozzle* The 
two-phase mixture drives the turbine rotor. The vapor is condensed, and the 
condensate is pumped back to the heater by a pump on the turbine shaft. 

The cycle is shown on a T-S diagram in Fig. 3. The state points are 
nunbered to correspond with Fig. 2. The two-phase nozzle expansion takes the 
fluid from a saturated liquid at 150°C (point 1) to a two-phase flow of 0.6 
quality at 40^C (point 2). The flow is decelerated in the rotor, condensed 
to point 4, and pumped back to nozzle inlet pressure at point 5. The liquid 
is then reheated by the source fluid to point 1. 

The feature of this cycle that makes it ideal for recovery of energy from 
a sensible-heat source is the straight-line heating of the working fluid from 
point 5 to point 1 over the full temperature range of the cycle. The heating 
curve of the working fluid matches the cooling curve of the source fluid. The 
temperature difference between the source fluid and the working fluid can be 
kept small at all points in the heat exchanger. By contrast, the heating 
curve of a Rankine cycle, in which the working fluid must be boiled to dry 
vapor, is a poor match to the source-fluid cooling curve. Large differences 
are necessary between the source-fluid temperature and working-fluid tempera- 
ture over portions of the heat exchanger, and it may not be possible to cool 
the source fluid all the way to ambient temperature; both effects reduce the 
available work. 

The appropriate efficiency for comparison of waste-heat cycles is the 
"utilization efficiency," defined as the ratio of power output to heat 
available. The heat available is the heat that could be extracted by cooling 
the source fluid all the way to the low-temperature limit set by ambient 
Lttupeiftiure or other external limitations. '*'f T_ is this low-temperature 
cooling limit, then the available heat is 








and the utilization efficiency is 


n 


u 


P 

”'=p‘^A - ■'C> 


( 2 ) 


where P is the power output, and m is the flow rate, c^ the specific heat, 
and the initial temperature of the source fluid# 

If, due to cycle limitations, the gas is cooled only part way to the 

cooling limit and leaves at a temperature higher than T^, then the heat 

input to the cycle is me (T T.) and the cycle efficiency, as 

pa a 

conventionally defined, is 


n (3) 

“p‘^A - ■'b> 

From Eqs. (2) and (3), it can be seen that the utilization efficiency 
is the cycle efficiency multiplied by the ratio of actual to available source- 
fluid temperature drop< 


n 



(4) 


In Fig. 4, the utilization efficiencies of two-phase and Rankine cycles 

are compared for the case in which the cooling limit is the same as the conden- 

. o 

sxng, or rejection, temperature ■ 38 C, and the turbine and pump 

efficiencies are 100 percent. The two-phase cycle is SO percent more effi- 
cient than the Rankine cycle for source fluid temperatures of 150 to 200^0. 

The utilization efficiency of the two-phase cycle is close to the limitip" 
efficiency W^/Q^ corresponding to the available work given by inte- 
gration of the Carnot efficiency between T and T . 

A K 

In an effort to better match the working fluid heating curve to the 
source fluid cooling curve, some geothermal studies have proposed supercrit- 
ical Rankine cycles. Figure 5(a) shows such a cycle using isobutane as 
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Fig* 4. Comparison of Rankine aud iwo-phase cycles 
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Fig. 5* Comparison of supercritical and two-phase cycles: 
(a) supercritical cycle (Isobutane); (b) two-phase 
cycle (Refrigerant 113) 
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Che working fluid* The uCiliseCion efficiency, with ■ 38°C and 

100-percenC turbine and pump efficiencies, is 0.157. A two-phase cycle using 
Refrigerant 113, Fig. 5(b), has a utilisation efficiency of 0.173, a 
10-percent gain. 

At higher source-fluid temperatures, water could be used as the working 
fluid. The main problem with water in a two-phase turbine is the large steam 
volune and high velocity at low exhaust pressures. Figure 6 shows how water 
could be used at high pressure by expanding only to 150^C in the two-phase 
cycle and using an organic bottoming cycle for the reaiainder of the tempera- 
ture drop. In this way a good temperature match between the working fluid and 
the source fluid could be maintained at a source-fluid temperature as high 
as 350°C. 

Another application of two-phase organic turbines is in refrigeration. 
Instead of throttling a refrigerant to produce cold liquid, the refrigerant 
could be expanded in a two-phase turbine. The work conventionally lost in the 
irreversible throttling process trould be recovered as shaft power. The shaft 
power could assist in driving the compressor. The electricity consumption for 
refrigeration would be reduced by about 10 percent. 

C. Wet-to-Dry Cycle 

If the initial temperature of the working fluid is sufficiently high and 
the saturated vapor line has a positive slope on a T-S diagram (a "drying" 
fluid), the working fluid could be expanded all the way to dry vapor. This 
phenomenon opens up the possibility of a two-phase cycle that has two-phase 
flow in the nozzle but not in the rotor. Such a cycle will be called a 
"wet-to-dry“ or WD cycle. 

Figure 7 is a T-S diagram of a WD cycle using toluene. The toluene is 
expanded from saturated liquid at 290^C (point 1) to saturated vapor at 
66^C (point 2). The vapor drives an impulse rotor and leaves the rotor 
slightly superheated at point 3. The vapor is condensed to point 4 and pumped 
back to the nozzle inlet pressure at point 5. There is again an ideal match 
to the cooling curve of the source fluid, but two-phase rotor flow is not 
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Fig. 6. Ste«ni topping cycle with organic-fluid bottoming cycle 
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Fig* 7* Wet-to-dry cycle 
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required. The ucilisetion efficiency for Che cycle of Fig. 7 ie 0.19 whereat 
the ucilisation efficiency of Che beat toluene Rankine cycle uaing the aane 
aource fluid, and having Che tame turbine efficiency, %K>uld be only 0.14. 

D. TWo-ComponenC Cycle 

Two-phase mixCurea in which Che vapor and liquid are different chemical 
compounds are called "two-component" mixtures. An example is steam mixed with 
an organic liquid or a molten salt. Such mixtures behave quite differently 
from the one-component mixtures discussed to this point. 

Figure 8 shows Che equipment arrangement for a two-phase, two-component 
engine, and Fig. 9 is a T-S diagram for the component Chat forms the vapor 
phase, in this case steam. The liquid phase is a molten salt such as lithium 
carbonate or sodium hydroxide (it is not known if these liquids would actually 
be compatible with steam). The ratio of liquid to vapor is 30, by mass. 

The steam and liquid mixture expands from the nosale inlet, point 1, to 
Che steam condenser pressure at point 2. Because of the large ratio of liquid 
to steam, the mixture expands almost isothermally and Che nossle exit velocity 
is much lower Chan it would be for steam alone. The two-phase mixture 
drives the turbine rotor. The steam is separated from Che liquid, cooled to 
saturation temperature in a regenerative heat exchanger (point 4) condensed to 
point 5, pumped back to the nossle inlet pressure at point 6, and heated to 
point 7 in the regenerative heat exchanger. The liquid leaving Che turbine is 
pumped to the heater by a pump on Che turbine shaft. The heated liquid is 
mixed with Che condensate (water) in the mixer, and the water is vaporised by 
direct-contact beat transfer. 

Possible applications of this cycle are in turbine engines for solar 
power or transportation. The advantage of the two-phase, two-component cycle 
for these applications is Che ability Co achieve a low shaft speed in a small 
turbine engine. A steam turbine of 150-kW shaft power, for example, using the 
temperatures given in Fig. 9, would have a speed of about 60,000 rpm. The 
two-phase turbine would operate at only about 10,000 rpm. There is also an 
efficiency advantage. At Che Fig. 9 temperatures, a steam Rankine cycle would 


15 









TEMPERATURE, •€ 







have an afficiancy of about 0<28, wharaaa tha two'-'phaac cycla would hava an 
afficiancy of about 0.37 (for tha aane turbine efficiency), becauae the ateam 
ia continuoualy reheated in the two-phaae noctle. In addition, the two-phaae 
cycle allowa control of turbine apeed becauae the mixture ratio can be varied 
to change the noxtle exit velocity, a capability unavailable in atean turbinea. 
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IV. NOZZLES 


A. Construction 

The nozzle converts the heat energy of the two-phase mixture to the 
kinetic energy of a jet. The gas is accelerated by pressure, and the liquid 
is partly accelerated by pressure, but mainly by drag. The gas moves faster 
than the liquid to produce the drag. The velocity difference, or *'slip", also 
causes breakup of the liquid into small drops. 

The nozzle design used in the turbine tests is shown in Figs. 10 and 11. 
The liquid enters the inlet fitting, flows around a baffle, and enters a 
hexagonal array of 61 tubes of 6.35-mm outside diameter and 4.57-mm inside 
diameter* The tubes are aimed toward the nozzle throat. The tubes are clamped 
together at their exits by a tube cage and are held 1.4 mm apart at the inlets 
by a tube sheet. The liquid enters the nozzle through the tubes at a velocity 
of about 4 m/s. 

The gas enters the mixer housing, flows through openings in the tube cage, 
enters the spaces between the liquid-injection tubes, and flows into the nozzle 
through the cusp-shaped gaps between the tube exits. The gas enters at about 
30 m/s. 

The multiple injection tubes and gas orifices provide a uniform 
distribution of flow rate and mixture ratio across the nozzle inlet. Where 
the inlet flow is saturated liquid, the multiple jets also provide an initial 
separation of the liquid streams so that vapor forms uniformly across the 
nozzle; a small amount of vapor is fed to the gas inlet to give a starting 
quality of 1 to 2 percent. 

The nozzle has a converging-diverging shape. The inlet section has a 
10-deg wall angle and the exit section a 2*5-deg angle. The throat section 
has a very gradual area change with distance to minimize slip* 


19 




2 









f 

The nozzle is designed for operation with Refrigerant 22 expanding from 
saturated liquid conditions at 20^C and 0.9 MPa to atmospheric pressure. 

The throat diameter is 13.1 mm. The exit is cut off at a 20-deg angle for 
turbine tests. The nozzle diameter at the upstream end of the exit ellipse is 
27.6 mm. The expansion area ratio to that point is 4.4, the value calculated 
to be required for expansion to atmospheric pressure. 

The nozzle is also used for water-and-nitrogen tests with expansion from 
2.0 MPa to atmospheric pressure. Figure 12 is a 1-v s flash photograph of 
the water-and-nitrogen jet leaving the nozzle. The velocity is 80 m/s. Faint 
striations from the liquid injection tubes are still visible in the jet. 

B. Nozzle Program 

Reference 10 presents the theory for two-phase nozzle flow. The theory 
is incorporated in two computer programs, one for single-component flow and 
the other for two-component flow. The programs, test cases, and operating 
instructions are available from JPL. 

To use the computer program, the nozzle inlet conditions are specified 
(pressure, temperature, flow rate, and velocity). The variation of pressure 
with distance and the fluid properties are entered in tabular form. The 
program calculates the liquid and gas flow rates, liquid and gas velocities, 
nozzle area, and other quantities at successive, small pressure steps. The 
results are printed at the exit pressure and any desired intermediate 
pressures. 

To match a given nozzle shape, rather than a given pressure profile, the 
specified pressure profile can be modified in two or three successive computer 
runs to make the calculated area variation agree with the desired nozzle 
contour. 

The liquid drop diameter is calculated from a Weber number criterion. 

The Weber number is proportional to the ratio between the pressure force 
tending to break up the drops and the surface tension force holding the drops 
together. The Weber number is defined as 
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o.(V - V )^D 
& R I 


where is the gas density, is the gas velocity, is he liquid 
velocity, D is the drop diameter, and a is the surface tension. 


Based on the drop breakup data reviewed in Ref* 10, the maximum stable 
drop diameter is the diameter for which the Weber number is 6* From Eq* (5), 
the corresponding maximum drop diameter is 


g g ^ 


In the nozzle program, the drop diameter D is initially set to a large 

value, and the diameter is reduced to D whenever D falls below D. 

* max max 

Typically, breakup takes place over a short distance immediately upstream of 
the throat. 


In the one-component nozzle program, the liquid and vapor are both assumed 
to have a temperature equal to the saturation temperature at the local pres- 
sure. In the two-component program, the temperature difference is calculated 
from drop heat transfer relations. 

Wall friction is calculated using nozzle friction relations for single- 
phase flow evaluated at the two-phase mixture density. 

C. Performance Definitions 

In a two-phase nozzle, the mixture ratio R is defined as the ratio of 


liquid mass flow rate m^ to gas flow range m . 
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The quality x is the ratio of gas flow rate to total flow rate 





( 8 ) 


where & is the total flow rate m. m . 
t t g 

The maximum velocity thermodynamically possible for expansion from the 
inlet conditions to the exit pressure is the isentropic velocity The 

mass weighted mean velocity of the actual flow with differing 
liquid and gas velocities is 


V - 


m V 
i I 


♦ mV 


m 


t 


(9) 


where is the liquid velocity and is the gas velocity. The mean 
velocity V is also equal to the jet momentum per unit mass. 

In this report, V will be used for the mean velocity with wall friction 
included, and will be used for the mean velocity without wall friction. 


The thrust of a two phase nozzle, by the definition of V, is 

F - m^V (10) 

Equation (10) is used to calculate the experimental value of the mean jet 
velocity V; the thrust and flow rates are measured, and the thrust is divided 
by the total flow rate to give V. 


* In Ref. 10, V and Vq are denoted and V, respectively. 
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The velocity coefficient of a nossle ia the ratio of mean velocity to 
iaentropic velocity 



( 11 ) 


The kinetic power in the jet is 


P. ■ ^(in V 
jet It 


♦ m V‘ ) 
g g 




The iaentropic jet power is 


P. « 
1 


I • .,2 


The nozzle efficiency is 


n 


n 



( 12 ) 


(13) 


(14) 


De LLL Steain-and-Water Data 

T. We Alger at LLL measured velocities for s tearo-and-vater nozzles, and 
the results are presented in Ref# 11. Nozzle 2 of Ref# 11 will be analyzed 
here. Nozzle 2 had a throat diameter of 6.4 mm, an exit diameter of 31#8 mm, 
a diverging length of 60.5 mm, and a divergence half-angle of 12 deg. The 
convergence half-angle was 45 deg and the throat had a sharp corner. In 
nozzle tests at an inlet pressure of 2.41 MPa, an inlet quality of 0.129, and 
an exit pressure of 30 kPa, the measured velocity coefficient was between 
0.90 and 0.95 (Fig. 10 of Ref. 11). 

In Ref. 12 Alger presents measurements of drop diameters for steam-and- 
water nozzles. He used a light-scattering technique. The nozzles had thin 
rectangular cross sections with side-wall contours similar to the nozzles of 
Ref. 11, but the inlet pressure was reduced to 1.0 MPa and the exit pressure 


l 
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to 20 kPa to provide greeter jet treneperency • The aeeaured mees-medien drop 
diaoeter (which is about the right diameter to use for nossle performance 
calculations) was found to be about 2.4 un (Ref. 12, p. 107). 

Using the measured pressure profile for Nozcle 2 from Fig. 7 of Ref. 11 
as input, Fig. 13 compares the data from Refs. 11 and 12 with predictions of 
the JPL nozcle program. The nozzle program, using Eq. (6) to calculate the 
drop diameter, predicts that the drops will break up to a 5.5-ub diameter. 

The highest measured velocity coefficient of 0.95 corresponds to a drop 
diameter of 2.8 urn, in good agreement with the measured drop diameter of 2.4 um 
from Ref. 12. Thus the drop diameter predicted by the nozcle program is 
perhaps as much as a factor of 2 too large. 

E. Water-and-Nitrogen Data 

Figure 14 compares the theoretical and experimental exit velocities for 
water-and-nitrogen flow in the nozzle of Fig. 11. The theoretical mean 
velocity is 4 percent below the measured velocities, again indicating 
overest iroation of drop size in the program. 

The comparison given in Table 1 is at the mixture ratio of 68 used in the 
turbine tests. The column headed "Fitted" shows the program results fitted to 
the measured data to give the best estimates of the quantities that could not 
be measured. 

The actual throat area is 9 percent higher than the throat area 
calculated by the nozzle program; this means that the measured flow rate is 9 
percent below theoretical for a given throat area. During expansion, the 
water has a slight reduction in flow rate to 3.603 kg/s because of 
evaporation, and the flow rate of the gas phase increases correspondingly. 

The theoretical velocity ratio of gas to liquid is 1.46. It is 
reasonable to assume that the velocity ratio in the experimental nozzle is the 
same. The "fitted" velocities that have that ratio and agree with the 
measured mean velocity V of 94.3 m/s are 93.7 m/s and 137.0 m/s. 

The fitted nozzle efficiency is 0.782, which is 8 percent higher than that 
predicted by the nozzle program. 
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Tabic 1. Coapariaon of theoretical and experiaental noxtle 
perforaance using water and nitrogen 




Value 


Itea 

Measured 

Nozcle 

program 

Fitted 

Inlet 




o 

Temperature, T^ C 

22 

22 


Pressure, Pp kPa 

2000 

2000 


Total flow rate, a^ , kg/s 

3.657 

3.657 


Liquid (water) flow rate, m^ , kg/s 

3.604 

3.604 


Gas (nitrogen) flow rste, a^, kg/s 

0.053 

0.053 


Mixture ratiog R 

68.0 

68.0 


Velocity, V^, m/s 

3.6 

3.6 


Throat 




Pressure, , kPs 


920 


Mean velocity, , m/s 


57.8 


A A 2 

Area, , mm 

135 

124 


Exit 




Pressure, i kPa 

98.6 

98.6 


Liquid temperature, 


21.3 


o 

Gas temperature, , C 


15.3 


Drop diameter, d^, urn 


130 


Area, A^ , mm^ 

597 

403 


Liquid flow rate, m^ , kg/s 


3.603 

3.603 

Gas flow rate, m , kg/s 
R 


3.054 

0.054 

Mixture ratio 


66.7 
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Table 1 (coned) 


Value 

lean 


Noaale 

progran 

Pitted 

Mean frec-eCrean velocity, Vq, ■/• 


95.1 


Mean velocity, V, ■/$ 

94.3 

90.5 

94.3 

Liquid velocity, V^, ■/• 


89.9 

93.7 

Cat velocity, V , ■/• 
S 


131.5 

137.0 

Jet power, kW 


15.03 

16.32 

leentropic velocity, , n/e 


106.8 


Isentropic power, P. , kW 


20.86 


Velocity coefficient, 

0.883 

0.847 

0.883 

Efficiency, r\ 

n 


0.721 

0.782 

Flow area ratio, R, 


8.9 
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Tht quantity that chiafly influancaa turbina rotor afficiancy ia tha 
ratio of gaa flow araa to liquid flow araa; tha largar tha ratio tha 
lowar tha rotor afficiancy. Tha valua of for thia noaala ia 8.9, a 
•■allar valua than in practical turbinaa, but usaful for rcaaarch purpoaaa in 
■aking liquid affacta pradoaiinata. 

F. Rafrigarant~22 Data 

Tha nocala of Fig. 11 waa tastad with aaturatad liquid Refrigarant 22. 

The raaulta are preacntad in Table 2. 

A anall amount of superheated Ratrigarant-'22 vapor waa fad into tha gas 
inlet to give a starting quality of 2 percant. For tha measured total flow 
rata, tha nossle program calculated a throat area 4 percant lass than actual. 
The measured exit velocity is 2 percent higher chan theoretical. The fitted 
noxxle efficiency shows a 14-percent improveawnC over Che water-and-niCrogen 
nossle efficiency of 0.78. The efficiency is higher because the higher 
molecular weight of Refrigerant-22 vapor gives lower slip velocity, and the 
lower surface tension of RefrigeranC-22 liquid gives smaller drop diameter. 

C. LLL Turbine Nossle 

The nossle used in the Lawrence Livannore Laboratory (LLL) turbine tests 
(Refs. 13 and 14) was analysed by Che nossle program, and the results are 
presented in Table 3. 

The theoretical nossle efficiency is 0.865. The thc*«retical velocity 
coefficient is 0.90. The measured velocity coefficient (Table 16 of Ref. 13) 
was 0.94 when the nossle was tested without the exit duct used in turbine 
tests and 0.91 with the exit duct. The velocity coefficient was only 0.83 
based on thrust measurements during turbine tests. Because of the uncertainty 
in the nossle performance, the theoretical performance values given in Table 3 
will be the ones used in analysing the turbine tests. 
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Table 2 (contd) 


Item 

Meaaured 

Value 

Nozsle 

program 

Fitted 

Mean free-stream velocity » m/a 


125.3 


Mean velocity, V, m/a 

126 

123.1 

126 

Liquid velocity, V^, m/s 


118.0 

121 

Gas velocity, , m/s 


136.6 

140 

Jet power, kW 


10.20 

10.72 

Isentropic velocity, V^, m/s 


134.4 


Isentropic power, , kW 


12.09 


Velocity coefficient, 

0.938 

0.916 

0.938 

Efficiency, 


0.844 

0.886 

Flow area ratio, R 


102 


• a 
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Table 3. Theoretical perfonaance of LLL turbine noccle 


Item 

Value 

Inlet 

Temperature, T^ , °C 

224.6 

Pressurep , kPa 

2528 

Total flow rate, m^ , kg/s 

0.596 

Liquid (water) flow rate» p kg/s 

0.512 

Gas (steam) flow ratep m p kg/s 

0.084 

Quality, 

0.141 

Velocityp p m/s 

15.0 

Throat 

Pressure p p kPa 

1787 

Mean velocityp p m/s 

131.9 

Aredp , mm^ 

72.5 

Quality, x^. 

0.174 

Exit 

Pressure p p^ , kPa 

13.2 

Temperaturep 

51.3 

Drop diameterp d^p um 

22 

Areap A^ p mm^ 

2615 

Liquid flow ratep m^ p kg/s 

0.397 

Gas flow rate, , kg/s 

0.199 

Quality, 

0.334 

Mean free-stream velocity, V^, m/s 

633.6 

Mean velocity, V, m/s 

627.0 

Liquid velocity, V|, m/s 

508.7 
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Table 3 (contd) 


Item 

Value 

Gas velocity, V^, m/s 

863.0 

Jet power, kW 

125.5 

Isentropic velocity, V., m/s 

697.7 

Isentropic power, P^, kW 

145.1 

Velocity coefficient, 

0.899 

Efficiency, 

0.865 

Flow area ratio, R 
a 

3300 


! 



The drop diameter of 22 um calculated for the LLL turbine nosile ia 
several times larger than that for the small nozzles of Refs. 11 and 12, 
because the turbine nozzle is longer and the pressure gradient is smaller, 
resulting in less shear to break up the drops. 

H. Wet-to-Dry Nozzle 

A toluene wet-to-dry (WD) nozzle was analyzed to determine the theoreti- 
cal effect of expanding to dry vapor. The calculated efficiency is very high. 
For toluene expanding from 289^C and 0.01 quality to 70^0, the exit quality 
is 1.0 and the efficiency is 0.98. 

Figure 13 shows how the nozzle diameter, quality, and phase velocities 
vary from inlet to exit for a toluene WD nozzle of 1.0-kg/s flow rate. The 
specified pressure profile is shown in the upper curve. The large nozzle 
length of 0.76 m is chosen to give a 2.5-deg divergence half-angle. The 
throat diameter is 9.8 mm, the exit diameter is 53 mm, and the expansion area 
ratio is 29. The area ratio is large because the quality at the throat is 
only 0.24 and most of the vapor remains to be formed. The vapor exit velocity 
(representing the entire mass flow at the exit) is 503 m/s. The isentropic 
velocity is 508 m/s. 
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SECTION V 
ROTORS 


A* Erosion 

One of the uncertainties in two-phase turbines is blade erosion. Even if 
performance problems are overcome, will satisfactory lifetime be possible? 
There should be no problem at the low velocities of organic working fluids, 
but there may be a problem with steam-and-water mixtures. Blade erosion 
measurements in steam turbines have been correlated by the following relation 
(Ref. 15): 


mass loss per unit area ^ _ V)^*^ (15) 

mass of liquid impinging per unit area c 

where V is the speed of the impinging liquid drops and V is a threshold 

c 

velocity below which little or no erosion occurs. 

The threshold velocity given in Ref. 15 for 12-percent chrome steel is 
120 m/s. The turbine blades in a two-stage, two-phase turbine travel at about 
two-thirds the speed of the jet. The jet velocity can, therefore, be 360 m/s 
without erosion, by this criterion. However, the velocity for steam-and-water 
mixtures is higher than this at many conditions of interest; the liquid 
velocity for the LLL turbine, for example (Table 3), is 509 m/s. Therefore, 
rapid erosion is a possibility in steam-and-water two-phase turbines. 

For organic working fluids, the velocity is much lower. Refrigerant 113 
expanding from 150^C to 46^C and 0.64 quality in a waste-heat cycle, for 
example, reaches only 158 m/s. 

Erosion was also of concern in the liquid-metal MHD project (Ref. 2). 

Erosion tests were made with a water-and-nitrogen jet at 137 m/s. A blunt-tipped 
aliininum cone placed in the jet eroded back 1.8 mm in 1100 hours, but 304 
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stainless steel eroded at only a tenth of that rate. These results indicate 
that erosion rates of reasonably hard alloys at relative speeds below about 
120 m/s should not be a problem. 

A WD-cycle turbine would, in principle, escape erosion entirely because 
only dry vapor enters the rotor. 

B. Liquid Path 

The key fact about the flow leaving the nozzle and entering the rotor in 
a two-phase turbine is that the drops are large, the gas density is low, and 
the drops travel in a straight line like bullets until they hit something. 

Figure 16 shows the results of trajectory calculations for water drops in 
a water-and-nitrogen jet at atmospheric pressure. The mixture travels through a 
90-degree bend. Small drops follow the gas with only a small outward drift, 
and only the drops starting near the outer wall impinge. Large drops travel 
straight into the wall. The fraction of entering liquid that impinges on the 
outer wall is plotted against drop diameter for various channel curvatures and 
flow velocities. Even for the largest channel radius (50 mm) and smallest 
velocity (75 m/s) there is complete impingement of drops larger than 13 iim 
in diameter; a 13-uro drop entering at the inner wall would just strike the 
outer wall at the end of the turn. Even 2.5-vjm drops would suffer 
50-percent impingement (those drops entering the outer half of the channel) at 
150- to 250-m/s velocity and a 25-mm channel radius. Thus, drop diameters can 
be only 1 to 2 pm for small impingement fractions at practical conditions. 

Such diameters are below the range produced by two-phase nozzles. 

Liquid flow behavior was studied photographically in the turbine rotor 
obtained from LLL, which will be referred to as Rotor 1. The blade shape is 
shown in Fig. 17. Figure 18 is a 1-ps flash photograph of the flow leaving 
the rotor in water-and-nitrogen tests. The total flow rate is 3.6 kg/s, the 
mixture ratio of water to nitrogen is 45, the jet velocity (mean velocity V) 
is 80 m/s, the rotor speed is 1540 rpro (the speed for maximum efficiency), the 
blade speed is 42 m/s (downward in the photograph), and the shaft output power 
is 8.1 kW. The drop diameter calculated by the nozzle program is 150 pm. 
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FRACTION OF DROPS IMPINGING 



Fig. 16. Effect of drop diameter on liquid impingement 


A1 
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Figure 18 ehowt that all of the water haa iapingad on the concave 
(upward-facing) aidaa of the bladaa and formed thin liquid ahaata. The ehaata 
are only about 0.5 ^m thick, 5 percent of the blade apacing. The leading edge 
of each ahaat aaMrgea froai a blade whan the blade haa traveled about half way 
down acroaa the noetla exit; at that poaition the firat liquid intercepted by 
the blade haa had tiaM to croaa the blade* Siaiilarly, the trailing edge of 
each liquid ahaat laavaa when the blade haa traveled about a half noxcla width 
beyond the end of the nosale. 

To the eye, the flow haa the aaaM appearance aa the jet leaving a 
noxala* The rotor haa owrely alotred and deflected the noacle jet* 

Figure 19 ia an axial view of the aaine flow* The firat ahaat of exit 
flow appec-a oppoaite the middle of the noaxle exit* About three leading 
edgea can be aeen before the emerging aheeta are hidden behind the foreground 
flow* The liquid leavea in a direction tangent to the rotor but deflected 
slightly outward from the original direction of noxtle flow* 

The vector velocities, both relative and absolute, of the liquid leaving 
the rotor can be measured from the photographs* The exit velocity relative to 
the rotor ia in the plane of the liquid aheeta, and is at the speed of the 
advancing leading edges. The absolute exit velocity is in the direction of 
the envelope of the sheets at the velocity of advance of the sheets in that 
direction. The relative exit velocity measured from Fig* 18 ia 27 m/a (66 
percent of the relative inlet velocity) at 43 deg above horixontal (the 
water does not follow the blade contour to the full 49-deg exit angle). 

The absolute exit velocity ia 27 m/a at 46 deg below horixontal. 

Two undesirable flow effects can be seen in Fig* 18. Soote of the water 
leaves the rotor on the noxxle side, and the trailing edges of the liquid 
sheets do not detach cleanly* These effects, which reduce efficiency, will be 
discussed in a later section. 

Figure 20 shows the exit flow at a lower liquid flow rate and higher 
velocity. The liquid flow rate is 1.9 kg/a, the mixture ratio ia 11, the jet 
velocity is 141 m/a, the rotor speed ia 2770 rpm, the blade speed is 76 m/s. 
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and the shaft power is 13.8 kW, The drop diameter from the nozzle program is 
75 urn. Figure 20 shows again that the water leaves as thin sheets, but some 
mist is also visible above the sheets. Evidently, some of the drops from the 
nozzle are small enough to remain entrained in the nitrogen flow. The mist 
leaves the rotor in a more nearly horizontal direction than the liquid, 
showing that the nitrogen and entrained small drops suffer less velocity loss 
in the rotor than the liquid, and are ejected at a higher speed. The water 
velocity loss is greater than that in Fig. 18, as evidenced by the steeper 
downward angle of the exit flow. The relative exit velocity of the water 
sheets, as measured from the photograph, is 36 m/ s (51 percent of the relative 
inlet velocity), and the absolute exit velocity is 56 m/s at 63 deg below 
horizontal . 


Figure 21 shows the exit flow at a still lower mixture ratio. The liquid 
flow rate is 0.54 kg/ s , the mixture ratio is 4, the jet velocity is 152 ni/s, 
the rotor speed is 3260 rpm, the blade speed is 89 m/s, and the shaft power is 
5.7 kW. The drop diameter from the nozzle program is 64 urn. More mist is 
visible than in the previous photograph. The relative exit velocity of the 
water sheets is 40 m/s (57 pe^'cent of the relative inlet velocity), and the 
absolute exit velocity is 66 m/s at 66 deg below horizontal. 

C. Rotor Program 


The picture that emerges from drop trajectory calculations and flew 
photographs is that the bulk of the liquid entering a two-phase turbine 
travels straight into the blades while the gas travels through as though the 
liquid were not present. Therefore, to calculate the performance of a two- 
phase turbine rotor, it should be sufficient to model the liquid flow as 
straight-in impingement, unaffected by the gas, followed by film flow along 
the blade surface. The gas flow can be modeled as conventional gas turbine 
flow, unaffected by the liquid. 


Such a model is derived in Appendix A. The liquid flow is 
detail, but the gas is merely assumed to exert a known fraction 
ideal gas-phase torque at any given speed. The ideal gas-phase 
blades at zero turbine bpeed is 2 m V , because the gas flow is 


modeled in 
n of the 
force on the 
idea 1 ly 
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reversed without velocity loss. The force ideally decreases with speed in 
proportion to (1 - where is the blade speed. Thus, the 

gas-phase force is 


F = 2nraV(l-V. /V) 
g g g g b g 


The gas-phase efficiency ^ peak value of at a blade 

speed of half the gas speed. Thus, the specified gas torque factor is 
also the peak gas-phase efficiency. 

The rotor flow model of Appendix A is incorporated in a computer program 

that was used for the theoretical calculations presented in this report. The 

inputs to the program are the nozzle exit dimensions; nozzle angle relative to 

the rotor, rotor speed; liquid and gas flow rates; liquid and gas velocities; 

blade chape; rotor diameter; liquid density and viscosity; gas torque fraction 

(peak gas-phase efficiency) n ; and (for windage torque calculation) gas 

g 

density and viscosity. 

Several different torques, powers, and efficiencies are calculated. 

First, the program calculates the total force F^ exerted on the blades by 
the two-phase flow. Multiplying this force by the rotor radius to the 
center of the blades gives the “blade torque'* L^. Multiplying this torque 
by the angular speed w gives the “blade power" P^. 


= L, fi) 
b 


The windage, or disc friction, torque L is calculated from the Mann 

w 

and Marston correlations (Ref. 16). Subtracting windage torque from the blade 
torque gives the theoretical “rotor torque" L . 


L = L. - L (1 

r b w 

product of rotor torque and rotor angular speed is the “rotor power* 

P = L fu ( ; 

r r 
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The ratio of the blade power to the jet power is the *'blade efficiency** 


n 


b- 


n 


b 



( 20 ) 


where is the jet power given by Eq. (12). This is the most useful 

efficiency for the comparison of experiment with theory. 


The ratio of rotor power to jet power is the **rotor efficiency** 


n 

r 



(21 ) 


The ratio of blade power to isentropic power will be called the 

**nozz le-blade efficiency*' n , 

^ nb 



1 


( 22 ) 


where n is the nozzle efficiency defined in Eq . (14). The nozzle-blade 
n y 'I 

efficiency is the efficiency that would be attained with a large number of 
nozzles so that windage loss is negligible relative to output power. 

The ratio of rotor power to isentropic power is the "turbine efficiency" 


n 


t 



(23) 


The turbine efficiency is also equal to the product of nozzle efficiency 
and rotor efficiency. 


n^=nn (24) 

t nr 
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The turbine efficiency is the only efficiency that counts in practical 
applications. The other efficiencies are useful for analysis purposes. 

D. Rotor Test Apparatus 

Figures 22 and 23 show the apparatus used in the two-phase turbine 
experiments. The test rotor is driven by a 60-kW electric motor through a 
variable-speed belt drive and a gear box. The gear box is an automobile 
transmission driven from the drive-shaft side. Rotor speeds can be adjusted 
from 700 rpm to 6000 rpm, but the highest used in flow tests was 3500 rpm. 

The rotor shaft can also be locked for zero-speed tests. The gear box is 
connected to the rotor through a rotating strain-gage torque transducer. 

The nozzle used in the turbine tests is the one shown in Fig. 11. 

Various throat diameters, exit diameters, and cut-off angles have been used. 

In some of the tests, as shown in Figs. 22 and 23, the nozzle was pivoted on 
bearings and restrained by a load cell for nozzle thrust measurement. 

The nozzle is fed through flexible hoses. In Fig. 22, the nozzle is 
connected for Refrigerant 22 tests. Liquid Refrigerant 22 is fed from a 
nitrogen-pressurized tank through a hand-operated valve to the liquid inlet of 
the nozzle. Vapor is fed from a heated Refrigerant 22 cylinder (not shown). 

For water-and-nitrogen tests, the liquid line is connected to water pumps, 
and the hand valve is used for controlling the water flow rate. In those 
tests, the nitrogen line is connected to the gas inlet of the nozzle, and the 
nitrogen flow rate is controlled by a pressure regulator. 

A water sump is located under the deck plates beneath the turbine. 

During turbine tests, the deck plates are removed and the flow leaving the 
rotor discharges into the sump. In the Refrigerant-22 tests, the liquid 
vaporizes on contact with the water and is exhausted from the building by 
blowers. 
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Fig. 23. Single-stage turbine 
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E. 


Measurement Procedure 


The main quantity of interest for comparison with theory is the blade 

torque L^. Blade torque is measured by setting the rotor speed at a desired 

value, with nozzle flow off, and recording the torque required to drive the 

rotor. Then the nozzle flow is turned on, the variable-speed drive is 

adjusted slightly to reduce the speed to the no-flow value, and the torque is 

recorded again. The torque without flow is the windage torque L^. The 

torque with flow, which is a torque exerted by the rotor on the drive system 

(causing the electric motor to act as a generator and feed power back into the 

power line) is the rotor torque L^. The sum of the two torques (measured in 

practice as the change in the recorder pen deflection between flow off and 

flow on) is the blade torque L, . 

b 


lade power is calculated from Eq . (17), = L^a> 


efficiency is calculated from Eq . (20), n. 


b jet 

for the jet power the "fitted" value obtained 
nozzle thrust, the calculated exit mixture ratio, and 
velocity ratio from the nozzle program, as illustrated 


, and the blade 
The value used 
from the measured 
the calculated exit 
in Table 1. 


F. Water-and-Ni t rogen Turbine Data 


Rotor 1 was tested with water and nitrogen at the nozzle conditions given 

in Table 1. The nozzle was mounted at an angle A of 20 deg on the right 

noz 

side of the rotor as drawn in Fig. 17. 

Figure 24 presents the measured blade torque as a function of rotor 
speed. The torque varies from 125 N-m at zero speed to 40 N-ni at 2250 rpm , 
and the torque follows a straight line. The speed at which the torque 
extrapolates to zero is J1300 rpm, corresponding to a blade speed of 91 m/s, 
which is close to the nozzle exit velocity. 


Figure 24 also presents the theoretical torque calculated by the rotor 
program (Appendix A). The gas torque factor n (equal to the peak gas-phase 
efficiency) is specified as 0.8; the true value of is probably in the 
range 0.6 to 0.9, but the nitrogen flow rate is so small compared with the 
water flow rate in these tests that the uncertainty in has negligible 
effect on the theoretical torque. 






The theoretical blade torque varies from 132 N-m at zero speed to 59 N-m 
at 1850 rpm. Beyond this speed, the rotor program calculates that the flow 
impinges on the backs of the blades. This effect is not modeled in the 
theory, and therefore the theoretical calculations are only carried to 1850 
rpm. 


Figure 25 compares the theoretical and experimental blade efficiencies 

The peak experimental efficiency is 0.631 at 1650 rpm. The peak 
theoretical efficiency is 0.^09. The agreement between theoretical and 
experimental efficiencies is sufficiently close to show that the theory 
accounts for the major lossec, h.it that there are additional losses amounting 
to about 10 percent. The possible additional losses are discussed in a later 
sect ion. 

Table 4 compares the theoretical and experimental performance at the 
peak-efficiency speed of 1650 rpm. The windage torque is six times greater 
than theoretical, but the theoretical windage torque from the Mann and Marston 
correlations is for a tightly-fitting housing, which gives lower windage 
torque than that for an unenclosed rotor. The theoretical values of relative 
inlet velocity and relative exit velocity (using the notation of 
Appendix A) show that there is a substantial liquid velocity loss even at the 
favorable operating conditions (low gas/liquid area ratio ) of this 
turbine. 

G. Ref r igerant-22 Turbine Data 

Rotor I was tested with Refrigerant 22 at the nozzle conditions of Table 
2. Because only a few runs of about 10 seconds duration could be made due to 
cost and tank-size limitations, the tests were made only at the theoretical 
peak-efficiency speed of 1880 rpm. 

The results for the run with the steadiest conditions and most reliable 
data are presented in Table 5. The measured blade efficiency is essentially 
the same as in the wa te r-and-n i t rogen tests. However, this efficiency is only 2 
percent less than the theoretical blade efficiency using an assumed gas torque 
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BLADE EFFICIENCY 7,., percent 



Fig. 25. Comparison of theoretical and experimental efficiencies 
of wa te r-and-n i t rogen turbine 
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Table 4. Comparison of theoretical and experimental turbine 
performance using water and nitrogen (Table I 
nozzle conditions) 


I tern 

Va 1 ue 

Rotor 

Measured 

program 

Nozzle angle, A , deg 
n 

20 

20 

Rotor radius, , mm 

263 

263 

Rotor speed, N, rpm 

1652 

1652 

Rotor torque, L^, N-m 

57.2 

66.5 

Windage torque, , N-m 

2.3 

0.4 

Blade torque, , N-m 

59.5 

66.9 

Rotor power, , kW 

9.90 

11.50 

Blade power, P. , kW 
b 

10.29 

11.57 

Rotor efficiency, 

0.607 

0.70A 

Blade efficiency, 

0.631 

0.709 

Nozzle-blade efficiency, 

0.494 

0.555 

Turbine efficiency, 

0.475 

0.551 

Relative inlet velocity, m/s 


53.3 

Relative exit velocity, , m/s 


29.6 


58 


Table 3. Comparison of theoretical and experimental turbine 
performance using Refrigerant 22 (Table 2 
nozzle conditions) 


I tern 

Value 

Rotor 

Measured 

program 

Nozzle angle, A^, deg 

15 

15 

Rotor radius, R^ , mm 

261 

261 

Rotor speed, N, rpm 

1880 

18^0 

Rotor torque, L^, N-m 

32.2 

34.0 

Windage torque, L N-m 

w 

1.7 

0.5 

Blade torque, , N-m 

33.9 

34. 5 

Rotor power, , kW 

6.34 

6.69 

Blade power, , kW 

6.67 

6. 79 

Rotor efficiency, 

0.591 

0.624 

Blade efficiency, 

0.622 

0.633 

Nozzle-blade efficiency, 

0.551 

0.561 

Turbine efficiency, 

0.524 

0. 553 

Relative inlet velocity, , m/ s 


72.6 

Relative exit velocity, , m/ s 


19.2 
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factor (gas phase efficiency) of 0,8. The uncertainty in the gas-phase 
efficiency introduces perhaps a 5-percent uncertainty in the theore^'cal 
torque because the liquid/gas mixture ratio is only 2.6. 

The measured turbine efficiency of 0.52 is the highest efficiency achieved 
with a two-phase turbine to dale. 

H. LLL S team-and-Water Turbine Data 

A two-phase turbine was tested by Lawrence Livermore Laboratory with a 
water-and-steam mixture at the nozzle conditions of Table 3. The results are 
reported in Refs. 13 and 14. 

The measured turbine efficiency of the LLL turbine is plotted in 

Fig. 30 of Ref. 13. Blade torque is recovered from those measurements as 

follows: The values are multiplied by the isentropic jet power ■ 

145.1 kW (Tab’e 3); the resulting rotor power values are divided by the 

rotor speed w at each data point to recover the rotor torque values L^; the 

experimental windage torque L^ at each data point is then calculated from 

the measured moment coefficients given in Fig. 29 of Ref. 13; adding L^ to 

L gives the measured blade torque L. ; multiplying L by rotot speed u) 
w ho 

gives the measured blade power P^. 

The windage torquo in the LLL ^ests was very large — 30 percent of the 
blade torque and six times the torque calculated from the Mann and Marston 
correlations for a closely fitting housing. Because it cannot be expected 
that the windage torque measured in motoring cests is exactly the same as the 
torque that should be added to rotor torque to obtain blade torque, there is 
an uncertainty of perhaps 5 to 10 percent in the recovered blade torque values 
for the LLL turbine. 

The theoretical blade torque of the L! L turbine is calculated using a gas 
torque factor n of 0.8, which introduces an uncertainty of several 
percent since the steam flow is 33 percent of the total. An even larger 
uncertainty is introduced by using the theoretical nozzle velocities from 
Table 3 as inputs to the rotor program since reliable nozzle measurements were 
not available. 
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Figure 26 compares the theoretical and experimental blade torques of the 
LLL turbine. The theoretical t >rque curve ends at 170C rpm, which is the 
speed at which liquid theoretically starts striking the backs of the blades. 
The torque results are similar to those obtained with the JPL water-and- 
nitrogen turbine, with the measured blade torques 10 to 20 percent below 
theoretical . 

Figure 27 compares the theoretical and experimental blade efficiencies. 
The highest mr-^sured efficiency is 0.47 at 1815 rpm, 10 percent below 
theoretical. With a slightly lower blade angle to avoid impingement on the 
back, the rotor program predicts that the blade efficiency would reach 0.62 at 
3200 rpm. 

Table 6 compares the theoretical and experimental performance in more 
detail at the highest test speed for which the water theoretically cleared the 
backs of the blades, 1575 rpm. This was also the speed at which the measured 
turbine efficiency n ^ reached its maximum value of 0.233; because of the 
large windage torque, the efficiency decreased again at higher speeds. 

A significant theoretical result shown in Table 6 is that the relative 
velocity of the liquid leaving the blades is only 17 m/s, 5 percent of the 
relative inlet velocity. This means that the torque could be calculated 
simply by assuming that the liquid stops on the blades and is thrown off at 
rotor speed. This approximation reduces the theoretical torque of the liquid 
by only A percent. The large velocity loss of the liquid is a consequence of 
the very high gas/liquid area ratio = 3300 (Table 3) which forces the 
liquid to spread over a large area in a very thin film. If the relative 
liquid velocity on the blade stayed at half the impinging velocity, the liquid 
film thickness would be only 22 pm. 

I. Staging Method 

Higher turbine efficiencies can be obtained by using two or more rotor 
stages. The approach is to reduce the velocity difference between the liquid 
and the blades and thereby reduce the energy dissipation due to friction. 
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Fig. 26. Comparison of theoretical and experimental torques 
of LLL steam-and-water turbine 







Table 6. Comparison of theoretical and experimental 
performance of LLL steam-and-water turbine 
(Table 3 nozzle conditions) 


I tern 

Measured 

Value 

Rotor 

program 

Nozzle angle, A^, deg 

20 

20 

Rotor radius, , mm 

973 

973 

Rotor speed, N, rpm 

1575 

1575 

Rotor torque, L^, N-m 

205 

329 

Windage torque, L , N-m 
w 

99 

17 

Blade torque, , N-m 

304 

346 

Rotor power, , kW 

33.8 

54.3 

Blade power, P^ , kW 

50.1 

57.1 

Rotor efticiency, 

0.269 

0.433 

Blade efficiency, n, 

b 

0.399 

0.455 

Nozzle-blade efficiency, 

0.345 

0.394 

Turbine efficiency, 

0.233 

0.374 

Relative inlet velocity, , m/ s 


362 

Relative exit velocity, , m/ s 


17 


6A 


Figure 28 illustrates the principle. The upper diagram shows a single- 

stage turbine blade in which the friction loss is so large that the flow is 

brought to a complete halt relative to the blade. The flow rate is unity, the 

jet velocity is 3, in arbitrary units, and the blade is traveling at the peak- 

efficiency speed of 1.5, in the same units. The flow is slowed from 3 to 1.5, 

exerting a force of 1.5. The blade power, equal to force times velocity, is 

2 

2.25. The jet power is 3 /2 or 4.5. Thus, the blade efficiency is 0.5. 

In the lower diagram, two stages are arranged in series so that the flow 
leaving the first stage enters the second. The flow is still brought to a 
halt relative to the blades in each stage. The first stage has a blade speed 
of 2. The flow is slowed from 3 to 2 in the first stage, exerting a force of 
1 and providing a power of 2. The first-stage exit flow enters the second 
stage at the velocity of 2 and is slowed to 1, exerting a force of 1 and 
providing a power of 1. The total power of the two stages is 3, giving an 
efficiency of 0. 66. 

With N stages, the efficiency is N/(N 1). Thus, even with large 
friction and complete velocity loss in each stage, the efficiency can, in 
principle, approach unity with a large number of stages. 

This same principle is used in drag-disc, or Tesla turbines, where the 
flow spirals inward between parallel discs, slowing due to the friction of the 
discs but always moving slightly faster than the disc speed. The efficiency 
can approach unity if the velocity difference between the fluid and the discs 
can be kept small throughout the inward spiral path. 

J. First-Stage Rotor Data 

To investigate staging experimentally, a second rotor was built, 
designated "Rotor 2," for use as the first stage with Rotor 1 as the second 
stage. Figure 29 shows the blade shape of Rotor 2. The blade angle is much 
shallower than in Rotor 1, because Kotor 2 runs at high speed, and the flow 
enters at a shallow angle. Before the final curved contour was machined into 
the blades, the rotor was tested with the flat blade shape indicated by the 
dashed line or the upper blade in Fig. 29. 
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I 


BLADE SPEED = 2 
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JT = 1 
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TOTAL POWER = 2+1=3 

EFFICIENCY = | 

A • D 3 


Fig. 28. Staging method: (a) single-stage turbine; (b) two-stage turbine 
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Figure 30 is a 1-ps flash photograph of the water-and-n i t rogen flow 
leaving Rotor 2 at 2200 rpm. The rotor has the flat blade contour in this 
photograph. The liquid leaves as thin sheets, as with Rotor 1. The downward 
angle of the exit flow is large because of the high rotor speed. Figure 31 
shows the same flow from the axial direction. The liquid can be seen reaching 
the blade exit on the near side after the blade has traveled about four blade 
spacings from the initial point of interception of the nozzle flow. 

The improvement in the flow pattern after the flat blades were remachined 
to the curved contour can be seen in Fig. 32. With the flat blades, the flow 
leaves at angles ranging from 57 deg above horizontal at the top of the jet to 
only 30 deg above horizontal at the bottom. With the curved blades, the flow 
leaves at angles between 60 and 45 deg. Evidently the curved blades provide 
needed centrifugal force during impingement to keep the liquid from bouncing 
back from the blade surface. 

A comparison between theoretical and experimental blade torques is made 
in Fig. 33 at the nozzle conditions of Table 1. Above 2300 rpm, the liquid 
starts impinging on the backs of the blades, and the theoretical curve ends at 
that point. The assumed gas torque factor is again 0.8, which may be 
optimistic for this rotor, but the uncertainty has little effect because of 
the small gas flow. 

The difference between theoretical and experimental torques is greater 
than that for Rotor 1. At 2200 rpm, the theoretical blade torque is 54 N-m 
and the measured blade torque is only 40 N-m, 25 percent less. The possible 
reasons will be discussed in a later section. The higher torque of the curved 
contour, especially at low speeds, is probably the result of the higher angle 
of the exit flow as shown in Fig. 32. 

K. Two-Stage Turbine Data 

The two rotors were assembled as a two-stage turbine. As shown in Fig. 

34, the first stage (Rotor 2) is connected to the variable-speed electric 
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Lg. 32. Comparison of exit flow from flat and curved blades 

(zero blade speed): (a) flat contour; (b) curved contour 
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drive. The second stage (Rotor 1) is connected to a water-brake dynamometer. 
The right side of Rotor I, as drawn in Fig. 17, faces the exit of Rotor 2. 

There is a gap of 0.5 mm between the rotors. 

To make a test run, the first-stage rotor is set at the desired speed. 

The flow is turned on and the dynamometer water flow is adjusted to give the 
desired second-stage speed. The torque of the first stage and the torque on 
the dynamometer are recorded, together with the rotor speeds. 

Windage torque is measured separately with the nozzle flow off. The 
windage torque exerted on the second stage by the first stage is sufficient to 
rotate the second stage for these measurements. 

The nozzle position can be adjusted radially and tangentially, and there 
is an optimum position for maximum total shaft power. Figure 35 shows the 
effect of moving the nozzle. The best radial position is R^ ■ 261 ram, at 
which position the nozzle exit ellipse is tangent to the outer wall of the 
blade passages, as shown in the sketch. The maximum second-stage torque occurs 
with a nozzle displacement 8 mm to the right of center (x ■ -8 mm), but the 
first-stage torque is low. The highest-power combination of first- and second- 
stage torque is with the nozzle located 25 mm left of center (x ■ 25 ran). 

In Fig. 36, the theoretical second-stage torque curve ends at 700 rpm 
where impingement on the backs of the second-stage blades begine. The 
theoretical torque is calculated using the theoretical exit velocity of the 
first-stage rotor. Because the measured torqu*^ is higher than theoretical, 
the absolute exit velocity from the first-stage rotor is evidently higher than 
theoretical. 

Figure 37 presents the efficiency of the two-stage turbine. The blade 
efficiency of 0.56 at zero second-stage speed represents the output of the 
first-stage rotor. The efficiency reaches a maximum of 0.70 at a second-stage 
speed of 750 rpm. The peak theoretical efficiency is 0.85. 

Table 7 compares the theoretical and experimental performance at a 
second-stage rotor speed of 678 rpm. The experimental windage torques are 
several times higher than theoretical but have the right sign (positive, or 
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Table 7. Comparison of theoretical and experimental 
two-stage turbine performance (Table 1 
nozzle conditions) 


Item 

Experimental 

First Second , 

Total 

stage stage 

1 

Theoretical 

First Second ^ ^ , 

Total 1 

stage stage j 

! 

Nozzle angle, , deg 

20 



20 



Rotor radius, , mm 

261 

261 


261 

261 


Rotor speed, N, rpm 

2200 

678 


2200 

678 


L 

Rotor torque, L^, N-m 

35.5 

33.6 


53.0 

21.4 


Windage torque, L , N-m 
w 

4.5 

-2.4 


0.5 

-0.2 


Blade torque, , N-m 

40.0 

31.2 


53.5 

21.2 


Rotor power, , kW 

8. 18 

2.39 

10.57 

12.21 

1.52 

13.73 

Blade power, P, , kW 
b 

9.22 

2.22 

11.44 

12.33 

1.51 

13.84 

Rotor efficiency, 



0.648 



0.841 

Blade efficiency, n, 

b 



0.701 



0.848 

Nozzle-blade efficiency, n 

no 



0.549 



0.664 

Turbine efficiency, 



0.507 



0.659 

Relative inlet velocity, , m/s 




42.5 

21.3 


Relative exit velocity, , m/ s 




31.3 

8.8 
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retarding, for the first stage and negative, or accelerating, for the second 
stage). The turbine efficiency is 0.51 experimentally and 0.66 theoretically. 

Although the measured turbine efficiency is not as high as that for the 
single-stage Refrigerant-22 turbine because of the lower nozzle efficiency, 
the blade and rotor efficiencies are the highest that have been achieved so 
fai with a two-phase turbine. 

L. Divergence and Stagnation Losses 

An additional source of loss not included in the theoretical torques and 
efficiencies is the radial spreading of the flow leaving the blades. To see 
the spreading better, a single blade was tested in the apparatus shown in 
Fig. 38. Sheet-metal blades similar in shape to those of Rotor 2 (flat 
contour) were held in a frame. A scoop collected the flow leaving all but the 
last blade. The flow passing behind the last blade continued on undisturbed 
except where it hit the support frame. 

Figure 38(b) shows the spreading of the liquid sheet leaving the test 
blade. The exit flow consisted of a central sheet of small divergence angle j 

and an outer fan on each side having about a 45-deg divergence angle. | 

The effect of the spreading exit flow is to reduce the circumferential 

component of the exit velocity by the cosine ot some average divergence 

angle. In the rotor program, the relative exit velocity is denoted V^. A 

corrected exit velocity can be obtained by specifying a divergence angle 

A , . such that 
d 1 v 

I 

V, = V- cos A,. (25) 

3 3 di V 

The theoretical zero-speed torque of Rotor 1 in the single-stage tests 

(Fig. 24) can be brought dcwn to the measured torque by using a divergence 

angle A of 30 deg. For Rotor 2 (Fig. 33), a divergence angle of 30 deg 
d 1 v 

lowers the theoretical zero-speed torque from 124 N-m to 119 N-m, whereas the 
measured zero-speed torque is only 108 N-m. However, as was seen from Fig. 

32, the angle of discharge of the liquid is less than the 60-deg exit angle of 
the blades. A line drawn through the estimated center of the exit flow in 
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Fig. 32(b) has an angle of only 40 deg above horizontal. If it is assumed 
that the liquid is leaving at 40 deg instead of 60 deg and is also spreading 
at an average divergence angle of 30 deg, then the theoretical zero''speed 
torque is reduced to 111 N-m, close to the measured value. 

A second loss not included in the theory is the stagnation and trapping 
of a portion of the liquid in the rotor. Observing a rotor in operation, it 
can be seen that some liquid continues to leave the rotor long after the rotor 
has passed through the jet. This effect can be seen in Fig. 39, which is a 
1-us flash photograph of the flow leaving Rotor 2 at 2250 rpm. Streamers of 
watei leave the outer wall of each blade passage after the rotor passes the 
jet. These streamers persist for about 90 deg of rotation. 

Further details of the main flow and the delayed stagnated flow can be 
seen in Fig. 40. An arm with four blades machined in the shape of those in 
Rotor 2 (flat blade contour) was rotated through the jet at 2200 rpm and 
photographed at seven successive positions- At the top position, the blades 
are just starting to enter the flow. In position 2, the liquid sheets are 
starting to emerge. In positions 3 and 4, the behavior of the flow that 
passes straight through the rotor can be seen most clearly: liquid is scooped 

out of the jet, moved to the near side of the rotor, and released to continue 
on in its original direction at reduced speed. 

In position 5, it can be seen that some liquid has remained with the 
rotor. In positions 6 and 7, this stagnated liquid continues to drain toward 
the outer walls of the blade passages and spin off, leaving the rotor at full 
rotor speed. 

The amount of liquid leaving the rotor in the main jet was measured by 
placing a scoop over the exit jet leaving Rotor 2. The outline of the scoop 
is shown in Fig. 41. The scoop collected 1.2 kg/ s of water when the nozzle 
water flow rate was 1.7 kg/s. The collected flow rate was constant for rotor 
speeds from 1250 rpm to 2200 rpm. Thus, 0.5 kg/s of water, or 30 percent of 
the flow, was missing from the main exit jet. Further measurements, obtained 
by moving the scoop, showed that, of this 0.5 kg/ s , 0.24 kg/s was thrown off 
the rotor in the first 30 deg of rotation beyond the jet, 0,07 kg/s was 
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Fig. 39. Exit jet and delayed drainage flow from Rotor 2 
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Fig. 40. Exit flow from bladed arm 
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Fig. 41. Rotor 2 exit scoop for flow measurement 
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thrown off in the next 90 deg, O.OS kg/s was thrown off the rotor on the 
nozzle side near the nozzle, and 0.14 kg/s was thrown off the rotor around the 
remaining circumference. 

If, in addition to the exit divergence effect, 25 percent of the flow is 
considered to be stagnated in the rotor and thrown off at rotor speed, the 
theoretical torque of Rotor 2 at 2200 rpm is decreased from 54 N-m to 40 N-m, 
in agreement with the data. Thus, the flow stagnation effect, combined with 
the exit flow divergence effect, is sufficient to account for the difference 
between theoretical and experimental torques. 

A major improvement in blade efficiency (from 0.70 to 0.85 for *^he 
two-stage water-and- nitrogen turbine) would be possible if the stagnation and 
divergence effects could be eliminated. This might be possible with more 
complex blade shapes. On the other hand, the undesired effects might reflect 
fundamental secondary-f low phenomena that cannot be avoided. Further study of 
these effects should be the next direction of research. 

M. Separator Turbine Concept 

The separator turbine (Refs. 3, 6, 7) is an alternative to the multistage 
impulse turbine for overcoming the friction problem of two-phase flow. The 
problem, in essence, is that the large flow area required by the gas phase 
requires large blade areas that are subjected to the large shear of the liquid 
phase. The multistage impulse turbine attacks this problem by reducing the 
relative velocity between the liquid and the blades, thus reducing the shear. 
The separator turbine attacks the problem by separating out the gas phase, 
thus reducing the required blade area. 

Figure 42 shows the basic arrangement of a separator turbine. The 
two-phase jet impinges on the inside wall of a rotary separator, causing the 
liquid to separate out as a liquid layer inside the separator drum. Ideally, 
the separator drum has the same speed as the jet. 

A second rotor, the liquid turbine, is pocitionea concentrically within 
the rotary separator. The liquid turbine carries one or wore U-tube scoops 
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with their inlets immersed in the liquid layer. The liquid enters the U-tubes 
and is discharged in the opposite direction. 

If the U-tubes travel at half the speed of the separator drum, the 
absolute velocity of the liquid leaving the U-tubes is zero. All of the 
kinetic energy of the liquid has been converted to shaft power of the liquid 
turbine. The gas leaving the rotary separator can, in principle, be sent 
through a gas turbine to recover the gas kinetic energy. 

In another variation of the separator turbine, the U-tubes can be 
replaced by diffusers that feed the liquid to a hollow shaft in the liquid 
turbine. If the liquid turbine is run at a speed less than half that of the 
separator drum, the liquid leaving the shaft will be pressurized. At the 
proper speed, the liquid can be returned to the two-phase nozzle, and the 
shaft power output will have been reduced by an amount equal to the pumping 
power for the liquid return. 

In principle, the only fluid friction in the separator turbine is in the 
U-tubes. Since the U-tubes carry only liquid, not two-phase flow, the wall 
area is small ^nd the friction loss is reduced accordingly. 

In practice, there is also impact loss of the two-phase jet entering the 
separator; this loss consists of radial impact on the separator drum and axial 
impact on the rear wall. There is also external drag on the U-tubes immersed 
in the liquid layer, and windage loss. Appendix B presents the theory of 
separator turbines, taking these losses into account. 

When the separator turbine theory is applied to specific cases, the 
theoretical efficiency of a separator turbine is no greater than t*' ^ theo- 
retical efficiency of a single-stage bladed turbine. Figure 43 compares the 
rotor efficiencies of separator turbines and bladed turbines for a particular 
s team-and-water mixture. The efficiency of the two-stage impulse turbine is 
plotted as a function of first-stage rotor speed. The efficiency of the 
separator turbine is plotted against liquid turbine speed; the separator speed 
is about twice as great. The efficiency of the liquid turbine is assumed to 
be very h igh : 0. 9. 
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At large diameters, the separator turbine is superior to the single-stage 
impulse turbine because the jet enters the separator with relatively small 
impact loss, and the assumed U-tube turbine efficiency is higher than the 
efficiency of the bladed turbine. However, at these large diameters, the 
windage loss is large. At the other end of the curve, at small diameters, the 
windage loss of the separator turbine is small, but the impact losses in the 
separator are large. The peak efficiency of the separator turbine is achieved 
at a liquid turbine speed of 14,000 rpm (separator speed of 23,000 rpm) • 

The basic concept of reducing friction by reducing liquid shear area is a 
valid appro/’:h, but it is not clear how to build a low-loss separator to 
implement the concept. 

N. Turbine Performance Examples 

The nozzle program, coupled with the rotor program, gives a reasonable 
prediction of two-phase turbine efficiency. The nozzle program gives nozzle 
efficiencies that are slightly too low, but the rotor program errs in the 
other direction. It is of interest to see what these programs predict in 
practical cases. 

A geothermal steam-water or ** total-flow** turbine has been one of the 
applications of interest. The theoretical performance of a 5-MW turbine at 
the flow conditions studied by LLL is presented in Table 8. The flow 
conditions are taken from Table 1 of Ref. 13. 

The nozzle program predicts a nozzle efficiency of 0.89 for nozzles of 
1.0-m length. For full admission, 15 nozzles are mounted around the circum- 
ference of a rotor of 0.7-m radius and 0. 1-m blade height. A single-stage 
turbine operates at 3600 rpm and has a rotor efficiency of 0.65; the turbine 
efficiency is 0.58. The turbine efficiency increases to 0.63 with two stages 
and to 0.66 with three stages. 
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Table 8. Theoretical performance of 5-MW steam-and -water turbine 
(LLL design co»'’*^ion8) 

Item Value 


Nozsle inlet 

Temperature, T^ , \ 

Pressure, p^, kPa 
Flow rate , , kg/s 

Quality, 

Nozzle exit 

Pressure, P 2 > kPa 
Temperature, T^ , 

Drop diameter, d^, wra 

3 

Gas density, o^, kg/m 
Area , , m^ 

Quality, x^ 

Isentropic velocity, , m/s 
Isentropic power, P^^ , MW 
Mean velocity, V, m/s 
Liquid velocity, , m/ s 
Gas velocity, , m/s 
Jet power, P. MW 
Nozzle efficiency, 

Flow area ratio, R 

a 

Rotors 

Radius, R^ , m 
Blade height, , m 
Gas-phase torque factor, n 


223.5 

2482 

30.0 
0.189 

12.4 

50.1 
28 

0.084 

0.144 

0.357 

742 

8.258 

681 

564 

891 

7.319 

0.886 

4150 

0.70 

0.10 

0.80 


i 

f 


I 
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Table 8 (contd) 


Item Value 


e-stage turbine 

Rotor speed, N, rpm 3600 

Relative inlet velocity, Vp m/s 329 

Rotor power, P^, MW 4.776 

Rotor efficiency, 0.653 

Turbine efficiency, 0.579 

Two-stage turbine 

First-stage rotor speed, , rpm 4770 

Second-stage rotor speed, , rpm 2320 

First-stage relative inlet velocity, , m/s 264 

Rotor power, P^, MW 5.209 

Rotor efficiency, 0.712 

Turbine efficiency, 0.631 


Three-stage turbine 


First-stage rotor speed, , rpm 5390 
Second-stage rotor speed, , rpm 3530 
Third-stage rotor speed, , rpm 1770 
First-stage relative inlet velocity, m/s 235 
Rotor power, P^, MW 5.426 
Rotor efficiency, 0.741 
Turbine efficiency, 0.657 
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The goal of the LLL program was a turbine efficiency of 0.70; this does 
not appear to be attainable with s tearo-and-vater mixtures. 

A higher efficiency can be achieved by using an organic working fluid and 
transferring the heat from the geothermal fluid in a binary cycle. Table 9 
presents the theoretical performance of a turbine using Refrigerant 113 (at 
conditions that correspond to a geothermal fluid of lower temperature and 
quality than those for the steam-and-water turbine of Table 8). 

The turbine efficiency is 0.67 with one stage, 0.71 with two stages, and 
0.72 with three stages. Part of the improvement in efficiency with 
Refrigerant 113 is due to the larger nozzle exit quality (0.64 instead of 
0.36), giving higher weight to the assumed gas-phase torque factor of 0.8. 
However, the liquid-phase blade efficiency is also improved. 

The results of the comparison shown in Fig. 44 are typical of what might 
be achieved ultimately with two-phase turbines. Single-stage turbines with 
unfavorable working fluids such as water and steam will have efficiencies 
limited to the mid- 50-percent range. Multistage turbines with more favorable 
working fluids might achieve efficiencies in the 70-percent range. Depending 
on the thermodynamic and practical gains offered by two-phase cycles in a 
particular case, such turbine efficiencies may be sufficient for useful 
applications. 
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Table 9. Theoretical performance of 5-MW Refrigerant-113 turbine 


Item 

Value 

Nozzle inlet 


o 

Temperature, , C 

149 

Pressure, p^, kPa 

1196 

Flow rate, , kg/s 

460 

Quality, 

0.01 

Nozzle exit 


Pressure, p^, kPa 

97 

Temperature, T^ , °C 

46 

Drop diameter, d^, um 

23 

3 

Gas density, o^, kg/m 

7.18 

Area, 

0.237 

Quality, 

0.641 

Isentrop'*.c velocity, , m/s 

174.2 

Isentropic power P^ , MW 

6.979 

Mean velocity, V, m/s 

168.6 

Liquid velocity, V^, m/ s 

158.0 

Gas velocity, , m/s 

174.5 

Jet power, P , MW 

6.551 

Nozzle efficiency, 

0.939 

Flow area ratio, R 

340 

Rotors 


Radius , R^ , m 

0.89 

Blade height , , m 

0.13 

Gas-phase torque factor, 

0.80 
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Table 9 (contd) 


I 

9 


Item Value 

Single-stage turbine 

Rotor speed, N, rpm 830 

Relative inlet velocity, V^, m/ s 89 

Rotor power, , MW A. 686 

Rotor efficiency, 0.715 

Turbine efficiency, 0.671 

Two-stage turbine 

First-stage rotor speed, , rpm 1090 

Second-stage rotor speed, , rpm 435 

First-stage relative inlet velocity, Vp m/ s 72 

Rotor power, P^, MW 4.94 7 

Rotor efficiency, 0.755 

Turbine efficiency, 0.709 

Three-stage turbine 

First-stage rotor speed, , rpm 1237 

Second-itage rotor speed, , rpm 638 

Third-stage rotor speed, , rpm 302 

First-stage relative inlet velocity, m/s 63 

Rotor power, P^, MW 5.055 

Rotor efficiency, 0.772 

Turbine efficiency, 0.725 
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TURBINE EFFICIENCY 77 percenr 


R 113 


WATER 


TOP CURVE 

REFRIGERANT 113 EXPANDING FROM 149“C, 
1200 kPa, 0.01 QUALITY TO 46<»C, 97 kPa, 
0.64 QUALITY. FLOW RATE 460 k^^s 

BOTTOM CURVE 


STEAM-WATER EXPANDING FROM 223“C, 
2480 kPa, 0.19 QUALITY TO 50“C, 12 kPa, 
0.36 QUALITY. FLOW RATE 30 kg/s 


NUMBER OF STAGES 


Fig. 44. Theoretical efficiences of 5-MW steam-and-water and 
Refrigerant-113 turbines 







VI, CONCLUSION 


Conclusions can be drawn in two areas: analysis capability for two*phase 

turbines and, to a lesser extent, practical usefulness of two-phase turbines. 

A. Analysis Capability 

The JPL two-phase nozzle program calculates velocities that are within 2 
to 4 percent of measurements, but the corresponding efficiency error is 4 to 8 
percent. The predicted velocities are lower than measured, probably due to 
overestimation of drop size. A more sophisticated drop size routine is needed. 
As an expedient, the critical Weber number could be reduced to a lower number, 
such as 3, that would fit existing data. 

The r^tor program derived in Appendix A gives an upper performance limit. 
Combined with the nozzle program, it probably shows the best that can be 
achieved with any given working fluid and operating condition. 

The measured rotor efficiencies have been 10 to 15 percent below the 
calculated values. The cause appears to be flow stagnation and divergence 
effects. It is not known if these losses can be reduced or if they are 
fundamental limitations. 

B. Application Prospects 

The efficiency of single-stage two-phase turbines will be only about 30 
percent, but this is sufficient for replacing the throttling steps in geo- 
thermal plants, refrigeration systems, and other systems that use irreversible 
flashing processes. A two-phase turbine would add power in these applications 
no matter how inefficient the turbine. However, the value of the added energy 
output must at least equal the cost of the turbine and associated equipment 
for the turbine to be economic. This requirement keeps many new energy 
devices from use, and is a stumbling block for two-phase turbines as well. 
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For applications where the two-phase turbine must compete with vapor 
turbines, there is a difficult efficiency requirement* The efficiency 
advantage of a two-phase cycle, say for waste-heat recovery or geothermal 
power where the two-phase cycle improves the matching to the heat source, is 
only about 15 percent. The efficiency of vapor turbines is about 80 percent. 
Therefore, a two-phase turbine must have an efficiency of more than 65 percent 
to give a net gain in cycle efficiency. This is about the upper limit of 
efficiency for two-phcse turbines, and is probably attainable only with 
organic working fluids. The picture that emerges is that two-phase turbines 
look promising in organic-fluid waste-heat or geothermal binary cycles where 
there is a significant thermodynamic advantage in using saturated liquid 
expansion and where the flow conditions are conducive to the best two-phase 
turbine efficiency. Two-phase turbines using steam and water expanding to the 
low pressures required in most applications such as geothermal may be ruled 
out by low efficiency and possibly erosion. 

Perhaps the most attractive prospect for early use of two-phase turbines 
is the WD cycle where two-phase flow is used only in the nozzle and dry vapor 
is used in the rotor. 
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PRECEDING PAGE BLANK NOT RLMED 

APPENDIX A. 

ROTOR MODEL 


A. Turbine Geometry 


The behavior of the flow in a two-phase turbine ^otor will be analysed 
for the geometry shown in Fig* A-I* 


A two-phase nossle of width W and height H (equal if nozsie is 

n n 

circular) delivers flow at angle A to rotor blades traveling at velocity 

nos ^ j 

V^* The jet contains liquid of flow rate m^ traveling at velocity 

V. and gas of flow rate m traveling at velocity V • 
t £ ^ SL 


The blades have an inlet section of shallow curvature where the flow 

impinges and an exit section of sharper curvature where the flow is turned. 

The inlet radius of curvature is R^ and the exit radius of curvature is 

R^. The blade ends with an optional straight extension of length L ^* 

2 r o ext 


A center plane can be dram thiough the point where the blade sn.rfaces 
pass through the axial direction. Th.« blade inlet angle is and the blade 
exit angle is ^ 2 % measured from the center plane. The inlet curvature 
ends at angle A^ past the center plane. 


The radius of the rotor at the center of the nozzle is R^ , and the 
outside radius of the rotor is R^. The blade spacing at radius R^ is S^. 


The flow parameter that has the dominating effect on liquid friction loss 
is the ratio of gas flow area to liquid flow area R^ , which is found as 
follows: The flow area occupied by liquid at the nozzle exit is 

where is the liquid density. The gas flow area is R^ titles the 
liquid area, and the total flow area is (1 R^) times the liquid area. 

Equating the product of (1 R^) and liquid flow area to the nozzle exit 

area (treating the nozzle as rectangular), the gas/liquid area ratio 

is given by 
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ANGLES POSITIVE 
TOWARD ARROWS 



i-ph<se turbine geoneCry 




1 


(A-l) 


R 


a 


WHO. 
n n t 



The liquid flow rate inpinging on each blade ia equal to the total 
liquid flow rate timea the fraction of the nosale area intercepted by one 
blade passage. The noszle has a projected width in the plane of the rotor 
equal to W /sin A • The fraction of the nozzle occupied by one blade 
passage is equal to the blade spacing divided by this projected nozzle 
width. Hence, the liquid flow rate per blade is 




m. S. sin A 
lb noz 


(A-2) 


B. Inlet Velocity Vectors 

Figure A-2 shows the velocity vectors for the liquid flow. Vectors are 
conveniently handled in the computer as complex numbers. A velocity V at angle 
A from the positive x axis forms a vector velocity V. If the magnitude V and 
angle A are given, the vector can be computed a^ V * POLAR(V, A), where POLAR 
is a function that constructs a complex number having real part V cos A and 
imaginary part V sin A. Conversely, given a complex number V, the magnitude 
and angle of V are given by V aBS(V) and A ■ PHASE(V), respectively, where 
ABS and PHASE are functions that operate appropriately on the complex number 
V. (Angles are in degrees in the computer program and will be written in 
degrees here. ) 

The two-phase jet is directed at angle 90 - A (degrees) below the 

noz 

X axis (A - 90 above the x-axis). The vector inlet velocity is thus 
noz 

V. - POLAR (V. , A - 90) (A-3) 

in t ’ fiCi. 

The vector blade velocity is 

V^ - POUR (V. , -90) (A-4) 

D D 


10 ! 





1 


The velocity of the incoming liquid relative to the rotor is the dif- 
ference between the jet and blade velocities: 



The magnitude of : .e relative velocity is 

= ABS (V^) 

and the angle of the relative velocity from vertical is 


(A-5) 


(A-6) 


A^^ = PHASE (V^) + 90 (A-7) 

noting that PHASE (V^) is negative. 

C. Impingement Geometry 


Figure A-3 shows the impingement geometry of a flow stream entering at 

distance Y above the blade inlet. The flow intersects the blade at angle 
s t 

A from the center plane. From Fig. A-3 it can be seen that angle A is 

X X 

given by 


j 

9 


sin A^ = sin A^ + (R^ cos A^ “ ^1 ^ ^st^ \l (A-8) 


Equation (A-8) can be rearranged for iterative solution as follows: 


A^ = sin ^ (C^ - tan A^^ cos A^) 


(A-9) 


where 

= sin A^ + (cos A^ - (A-10) 

The impingment angle of the streamline crossing at Y is then given by 

s c 

e = 90 - A„, - A (A-11) 

VI X 
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Impact Behavior 


D. 


When a single liquid jet strikes a surface, both momentum and energy are 
conserved. The jet splits into two streams that travel away from the point of 
impact in opposite directions at full velocity (conserving energy) and at 
differing flow rates (conserving momentum). 

When multiple jets or drops strike a surface, only momentum is conserved, 
because adjacent streams collide with each other and lose kinetic energy 
through inelastic impact. This effect is familiar from the household water 
faucets that use screens to divide the flow into multiple streams that dis- 
sipate their energy on impact and thus do not splash back. 

For two-phase flow striking a turbine blade, the process for the liquid 
phase can be described as shown in Fig. A-4. The liquid drops strike the 
blade at angle 0 and velocity Vy The row of drops striking the inlet 
edge of the blade (idealized as a sheet of liquid) splits into a forward flow 
and a smaller back flow, both flows traveling at velocity V^. The second row of 
drops also splits into a forward flow and a back flow, but the back flow 
collides with the larger forward flow from the first row of drops and is swept 
back into the forward direction. Each succeeding row of drops has its back 
flow swept into the forward direction, and the only back flow from the entire 
impinging stream is that from the first row of drops. 

For a large array of jets, the back flow thus carries a negligible 
fraction jf the total momentum, and the stream leaving the impingement zone 
has a momentum m that is essentially equal to the entire forward 
conponent of momentum of the incoming jet m cos 6. Thus, the equation 
for liquid velocity recovery in an impinging two-phase jet, considering impact 
only, is 


V 2 * cos 0 


(A-12) 
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E. Combined Impact and Friction Loss 

Figure A-5 shows a short section of blade having width H (that dimension 

n 

perpendicular to the plane of the figure) and area dA on which the flow ij 
impinging at angle e. The liquid flow rate impinging on the blade section is 
dm. The flow rate previously collected is , entering at velocity Vq, 

The flow rate in the film at the midpoint of the blade section is 

in * DU ♦ ^ dm. If the velocity of the film at the midpoint is V , then 
m U m 

the film thickness is 


t 

m 


m 

m 


PflH V 
t n m 


(A-13) 


It can be assumed that the wall friction is the same as for flow in a 

rectangular channel of height 2t and width H , because the velocity 

m n 

gradient at the film surface is zero as it would be at the center of a 
rectangular channel. The hydraulic diameter of a rectangular channel is four 
times the area divided by the wetted perimeter. Thus the hydraulic diameter 
of the film is 


4 (2t H ) 
m n 

2 (2t + H ) 

m n 


(A-14) 


Since t is much smaller than H. , the hydraulic diameter is 
m b’ 

approximately 


* 4t 
h m 


(A-15) 


The Reynolds number for pipe flow is 


Re 




(A-16) 


where is the viscosity of the liquid. 
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Substituting the hydraulic diameter from Eq. (A-I5) and the film 
thickness from Eq. (A-13)| the Reynolds number at the midpoint is 


Re 

m 



(A-17) 


For very low Reynolds numbers the flow is laminar and the friction 
coefficient is given by 


C 


f 


11 

Re 


(A-18) 


For turbulent flow a convenient expression for friction coefficient, 
valid over a wide range of Reynolds numbers, is the Von Karman equation: 

C • (A-19) 

(A log^Q (2 Rev/C^) - 1.6)*^ 

The Reynolds number at which Eqs. (A-18) and (A-19) give the same is 
1034.8, and this Reynolds number is used as the dividing line between the two 
equations. 

The retarding force on the liquid film in the blade section of Fig. (A-5) 
is the product of dynamic pressure, friction coefficient, and wall area. Any 
mean velocity V between inlet velocity and exit velocity V can be 
used for calculating dynamic pressure since the two velocities approach each 
other in the limit of small section lengths; exit velocity is used for 
V here. The friction retarding force can then be written 

F - *5 o.V^ dA (A-20) 

I e rm 

where Cf^ is the friction coefficient at the middle of the section. 
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The flow area of the two-phase stream impinging on the blade section is 

dA sin e. For unit liquid flow area the gas flow area is and the total 

flow area is 1 ^ R , where R is the gas/liquid area ratio. Hence, the 
a a 

liquid flow rate impinging on the blade section is 

0 | V. dA sin 6 

<14 . -I 4 -P 5 (A- 21 ) 


Solving Eq. (A-21) for area dA and substituting into Eq.(A-20), the 
retarding force on the liquid film in the blade section is 


V (1 + R ) dm 
e fm a 

2 sin 0 


(A-22) 


The momentum dM added to the film by the impinging flow dm, according to 
Eq.(A-12), is the momentum of the impinging stream times cos 0 , or 

dM * dm cos 9 (A-23) 

The momentum of the previously-collected flow entering the section along 
the surface is 


^0 * 


(A-24) 


The momentum of the flow leaving the section is 


m V 
e e 


(A-25) 


The exit momentum is equal to the sum of the inlet momentum and the 
momentun added by the impinging flow, less the friction force: 


♦ dM 


(A-26) 


Substituting the friction force F from Eq.(A-22) and the momentums from 
Eqs. (A-23, A-24, and A-25), the momentum equation becomes 
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cos 9 ■ 0 (A-27) 



(I ♦ R^) 

2 sin 9 






" 0 % 

dm cos 6 


Using the quadratic equation to solve for the ratio the result 

is 


V . ^ ra 

e sin 6 _e 

V “ C, ( I ♦ R ) dm 

1 fro a 





2 C 


iii.l_l.lii/ ^ (i 

tan e \^j \ 


" 0^0 

dm cos 6 


)- 


(A-28) 


This equation gives the liquid film velocity leaving a blade section 

on which two-phase flow or liquid flow rate din is impinging at angle 6 with 

liquid velocity V. and gas/liquid flow area ratio R , when the incoming 
1 a 

film flow rate is m^ at velocity and the outgoing liquid flow rate is 

m ■ m^ ♦ dm. 
e 0 


Equation (A-28) reduces to the equation used in past work (Reference 2, 

Figure 6, substituting Ajk * sin 6 and r ■ R there) for the 

z s V a 

velocity leaving a flat plate with no initial surface flow. In that case 
ro^ * 0, dm * , and the exit velocity is given by 


V 

e 



sin 6 

C- (1 + R ) 
f m a 



2 C. (1 ♦ R ) 
fm a_ 

tan e 


1 


(A-29) 


F. Film Friction Loss 

Following the impingement zone there is additional blade surface to com- 
plete the turning of the flow. The Reynolds number for this **f i Im-f low** zone, 
from Eq . (A-1 7) , is 


Re 


f 


4 ro. 


H u, 
n t 


(A-30) 
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where is the blade flow. The friction coefficient for this Reynolds num- 
ber as designated 

Figr A-6 shows the definition of various surface lengths along the 
blade. The length to the end of the inlet radius of curvature is 

h “ ^ ^3^ (A-31) 

The length of the impingement zone is 

L. - R, (A, - A ) (A-32) 

imp 1 1 X 

where A^ is the angle given by Eq . (A-9) for the last streamline of impinging 
flow (the flow envelope). 

The total surface length is 

■-.urf ■ h * “2<*2 - *3> * ‘•.xt <*-”> 

The length of the film-flow zone is 




surf 


L. 

imp 


(A-34) 


As the film decelerates, the friction force in any length dx is equal to 
the momentum decrease in that length. The riomentum change is 


“b dV - - *5 O^V 


(A-35) 


Integrating Eq. (A-35) over the length velocity decreases 

from V 2 to according to 


V 


3 


1 

''2 ' 2*1, 


(A-36) 


112 



c 


Exit Velocity Vectors 


The exit velocity vectors are shown in Fig. A-2. The vector exit velo- 
city of the liquid relative to the blade is 


V^ - POLAR (Vj, A^) 


(A-37) 


Adding the blade speed, the vector absolute exit velocity is 


V . ■ V- ♦ V. 
out 3 b 


(A-38) 


The magnitude of the absolute exit velocity is 


V - ABS (V ) 
out out 


(A-39) 


and the angle of the absolute exit flow from vertical is 


A - 90 ♦ PHASE (V ) (A-40) 

out out 

noting that PHASE is negative. 

For multistage turbines the exit velocity V ^ becomes the inlet velo- 

out 

city V. to the next stage, 
in 

H. Windage 


Mann and Marston (Ref. 16) present correlations for predicting the wind- 
age, or disc friction torque, on a bladed rotor. 

The Reynolds number of the flow on a disc of radius R^ rotating at 
angular velocity w is defined as 


Re 


„2 

0 w R 
g w 


where o is the gas density and u is the gas viscosity, 
g g 


(A-41) 
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If the rotor is one of several stages, then separate Reynolds numbers Rcj 
and Re^ are calculated for the upstream and downstream sides of each rotor, 
using the angular velocities relative to the adjacent rotors or to the 
atmosphere. 

The windage torque depends on the thickness of the rotor, which is equal 
to the blade chord. From Fig. A-3 it can be seen that the blade chord is 
given oy 


C - Rj(sin Aj + sin A^) + R 2 (sin A 2 “ sin A^) + cos A 2 (A-42) 

The aspect ratio of the blades is AR ■ where is the blade 

height (Fig. A-1). Table 3 of Ref. 16 gives moment coefficients C as a 

6 

function of aspect ratio for a Reynolds number of 10 and a chord/radius 

ratio C/R of 0.14. These coefficients are denoted C*. The table can 
w m 

be fitted by 


C„* - 0.0067 ♦ 0.007 ArO*^ 


(A-43) 


The exponent y of the Reynolds number dependence (C proportional to 
**y 

Re ) is given in the same table, and the values can be fitted by 


r - 0.21 - 0.086 AR 


(A-44) 


up to A?i ■ 0 . 7 , beyond which the ext rapolat ion is uncertain and y is held 
constant at 0. 1 3. 


The moment coefficient C for a disc of finite thickness is related to 

m 

the moment coefficient fot a disc of zero thickness by Eq . (3) of Ref. 
16: 


C 

m 


C (1 ♦ 2.3 

”0 V 



(A-4 3) 
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Using the exponent y to correct fron Re ■ «itid factoring out the 

thickness/radius ratio C/I^ ■ 0.14 for which Table 3 applies, the 
sero-thickness moment coefficient as a function of Reynolds number is 


C *(10^/Re)^ 
m 

m^ “ 1 ♦ (2.3)(0.14) 


(A-46) 


Equation (A-43) can then be applied to find the moment coefficient for 

the actual thickness/radius ratio C/R • 

w 

The windage torque M on one side of the disc is found from Eq. (2) of Ref. 

16: 

M - 0.25 C 0 (A-47) 

mg w 

For a single-stage turbine, n is the rotor speed and M is doubled to give 
the total windage torque L^. For multiple stages the relative w values on 
each side of the disc are calculated and used with the proper signs to give 
the windage torques and on the inlet and exit sides, respectively 
(positive for a retarding torque and negative for an accelerating torque). 

I. Torque, Power, and Efficiency 


The force exerted by the liquid on the blades is equal to the change in 
liquid momentum in the direction of blade motion. For the angles as defined 
in Fig. A-2, 


F ■ m (V. cos A. “ V ^ cos A 
t t in in out out 


(A-48) 


The torque of the liquid phase is the product of the liquid force and the 
rotor radius: 


L 


I 




(A-49) 


• 


The gas phase acts independently of the liquid phase, and the gas-phase 
torque can be calculated by conventional methods. The present program uses 
the approximation that the gas-phase torque varies linearly from maximum at 
zero speed to zero at synchronous speed (equal blade and gas speed), with a 
specified fraction of ideal torque denoted by n . The factor n is 
also the gas-phase efficiency at half synchronous speed. This assumption 
gives the following gas-phase torque equation 


2 n m V R^(l - VjV ) 
g g g b b g 


(A-50) 


The total blade torque is the sum of the liquid and gas torques, 
^ ^g * rotor torque is the blade torque less the 


windage torque, L 


H, - 


The blade output power is 




(A-51) 


and the net rotor power is 


P = L (D 
r r 


(A-52) 


These powers are summed over each stage for a multistage turbine. 


The power in the incoming jet is 


P. ^ + m V 

jet I I g g 


(A-53) 


The blade efficiency is the ratio of blade power to jet power: 


^jet 


(A-54) 


The rotor efficiency is the ratio of rotor power to jet power: 


r P. 


(A-55) 


jet 
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Tlie turbine efficiency is 


n 


t 


(A-56) 


where is the nozzle efficiency. 

J. Divergence and Stagnation Losses 


Equations (A-1) through (A-56) define what is meant by "theoretical** in 
this report. The losses modeled are impact and friction for the liquid , the 
specified torque factor n for the gas, and the Ref. 16 moment coefficients 
for windage. 


Two additional losses can be calculated in the program. The effect on 
torque o^ the spreading of the liquid leaving the blades is calculated by 
specifying a mean divergence angle for the exit liquid flow. The 

relative exit velocity is then corrected to 

V. * = cos A^. (A-57) 

j 3 div 

A fraction of the liquid can be specified as stagnating in the 

rotor and leaving at rotor speed. The force exerted by this liquid is 


F 

I 


stag 


m. f ^ 

I stag 


(V. cos A. - V^) 
in in b 


(A-58) 


Eq. 


The force of the main flow is (1 - f ) times the force given by 

stag 

(A-48). The total liquid force is that force plus . The result is 

stag 
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F. - m (V. cos A. - (1 - f ^ ) V cos A - f , V. ) (A-59) 

t t in in stag out out stag b 


The inlet velocity to the next stage is assumed to be the weighted 
average of the main-flow and stagnated-f low velocities: 


in 




(1-f^ )V^^f^ V. 

stag stag b^ 


(A-60) 


K. Optimization 

The program can be used in an optimization search mode. In this mode the 
program varies the rotor speeds and (optionally) the blade inlet angles for 
each stage until the maximum rotor efficiency is found. 

L. Program Listing 

A listing of the computer program follows. The subroutine ZPMIN 
referenced is a library routine for minimization. Decks of the rotor program 
including ZPMIN are available from JPL. 

The nomenclature used in the program is defined in Table A-1 . 
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TWO-PHASE ROTOR PROGRAM 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 


C 


C 

C 


• D. G. ELLIOTT, MARCH 31 » 19d0 • • 

• LAST REVISION NOV. 30, 1981 • • 

• • • • • ASCII FORTRAN •••#•• 
EQUATION NUMBERS REFER TO APPENDIX A 


PARAMETER MXSTG=2, MXSTCP:20 § MAX STAGES AND IMPINGQtENT STEPS 
PARAMETER N=MXSTG, N1=2*N, N2:N1*(NU3), N3rMXSTEP-«-1 

IMPLICIT REAL(L,M) 

REAL A1(N),A2(N),A3(N),R1(N),R2(N),RPM(N),LEXT(N) g INPUT ARRAYS 
REAL AIN(N),AOUT(N),AV1(N),CFF(N),FL(N),L1(N),LB(N) 

REAL LG(N) ,LIMP(N) ,LL(N) ,LSURF(N) ,LR(N) ,LW(N) ,REF(N) 

REAL VI (N) ,V2(N) ,V3(N) ,VB(N) ,VIN(N) ,VOUT(N) ,WW(N) 

REAL THETA(N3) ,X(N1 ) ,XBEST(N1 ) ,WA(N2) ,EPSX( 10) 

COMPLEX C V 1 ( N ) , CV3 ( N ) , C VB( N ) , CVIN ( N ) , C VOU T ( N ) , POL AR 

NAMELIST /IN/HN,HN,HB,AN0Z,RPM,ML,VL,MG,VG,R1 ,R2,LEXT,A1 ,A2,A3, 
& SB , RB , RW , RHOL , RHOG , VISCL , VISCG , NSTG , NSTEP , 

& ADI V , NR 1 , ETAG , ETAN . FST AG , NINPT , NDET , NO PT , NMO N 

NAMELIST /OUT1/LB,PB,ETAB,RPM 
NAMELIST /0UT2/LR,PR,ETAR,ETAT 
NAMELIST /OUT3/LL,LG,LW 

NAMELIST /OUT4/A1 ,VB, VIN,AIN,V1 ,AV1 ,V2,L1 ,LIMP,LSURF,V3 , 

& AOOT,VOUT,WW,RAP1 ,MB,REF,CFF,TFILM 

NAMELIST /OUT5/CV1 ,CV3,CVB,CVIN,CVOUT 
NAMELIST / 00 T6/M0 , MM , ME , REM , CFM , TERM2 , VO , VE 


••••• FUNCTIONS 


XSIN(A)=SIN(PX*A) € SINE OF ANGLE 'A* IN DEGREES 
XCOS(A)=COS(PX*A) § COSINE OF ANGLE 'A* IN DEGREES 
XTAN(A)=TAN(PX»A) € TANGENT OF ANGLE 'A' IN DEGREES 
PHASE(C)=180/PI»ATAN2(AIMAG(C), 

&SIGN(MAX(ABS(REAL(C)),1E-38),REAL(C))) wangle of C,DEG CCW from X-AXIS 
POLAR(AMP,ANG)=CMPLX(AMP»XCOS(ANG) ,AMP»XSIN(ANG) ) § VECTOR(AMPL, ANGLE) 

••••• CONSTANTS AND DEFAULT INPUTS ••••• 

DATA EPSF/1E-4/,EPSX/10»1E9/,ITMAX/100/ € CONSTANTS FOR SEARCH ROUTINE 

DATA A1/N»20./,A2/N»60./,NSTEP/20/ ^DEFAULT INPUTS 

PI=4.0«ATAN(1.0) 

PX=PI/180 
PXX= 180/PI 


120 


c 

c •§••• input ••••• 

c 

10 PRINT •,» ENTER* 

READ (5.1N,ERR:10,ENDs100) 

IF (NSTEP.GT.MXSTEP) PRINT NSTEP TOO LARGE* 

IF (NSTEP.GT.MXSTEP) GO TO 10 

C DEFINE FILE 6(APRINT,,72) « TO PRINT NAMELIST ON TERMINAL 
C 

IF (NINPT.EQ.1) WRITE (6, IN) 

C 

C ••••• INITIALIZATION ••••• 

C 

RAP1=WN«HN«RH0L»VL/ML # GAS/LIQ AREA RATIO +1(1) 
MB=ML«SB*XSIN(ANOZ)/WN « LIQUID FLOW PER BLADE (2) 

DM=MB/NSTEP § LIQUID FLOW PER BLADE SECTION 

CVIN(1)=POLAR(VL,ANOZ-90) § VCTR LIQUID VELOCITY FROM NOZZLE (3) 
PJET=0.5*(ML»VL»»2 ♦ MG»VG«»2) 9 JET POWER (53) 

C 

DO 20 K=1.NSTG € SET RPM GUESSES IF SEARCH USED AND SET VB*S 
IF (NOPT.GT.O) RPM(K)=30/(PI*RB)*VL«(NSTG+1-K)/(NSTG+1.) 

20 VB(K)=RPM(K)»PI»RB/30 9 INITIAL GUESS AT BLADE SPEEDS 

C 

IF (NOPT.EQ.O) GO TO 40 9 NO OPTIMIZATION 
C 

C OPTIONAL OPTIMIZATION OF RPM(S) AND BLADE INLET ANGLE(S) ••••• 

C 

DO 25 Ksl.NSTG 9 SEARCH INITIALIZATION 
X(K)=VB(K) 9 FIRST NSTG ELEMENTS IN X-VECTOR 
25 X(NSTG+K)=A1(K) 9 SECOND NSTG ELEMENTS IN X-VECTOR 

C 

KG0=0 

NX=NOPT»NSTG 9 NUMBER OF VARIABLES IN X-VECTOR 
30 CALL ZPMIN(NX,X,y,EPSX,EPSF,nMAX,WA,KGO,XBEST,YBEST) 9 RESET X 
C 

IF (KG0.GT.1) GO TO 90 9 CONVERGED OR ERROR 
C 

DO 35 Jsl.NSTG 9 COPY NEW X-VECTOR TO VARIABLES 
VB(J)=X(J) 9 NEW BLADE VELOCITIES 
35 A1(J)=X(NSTG+J) 9 NEW BLADE ANGLES 

C 

C #§••• STAGE CALCULATIONS ••••• 

C 

40 DO 70 K=1,NSTG 

VIN(K)=ABS(CVIN(K)) 9 ABSOLUTE INLET VELOCITY 
AIN(K)=PHASE(CVIN(K)) +90 9 ANGLE OF ABSOLUTE INLET VELOCITY 
CVB(K)=P0LAR(VB(K),-90) 9 VECTOR BLADE VELOCITY (4) 
WW(K)sVB(K)/RB 9 ANGULAR SPEED 
RPM(K)=30»WW(K)/PI 9 RPM 
C 

c §§#•• blade inlet conditions ••••• 

c 
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CV1(K)sCVIN(lO-CVB(K) « VECTOR EELATIVE INLET VELOCITY (5) 
V1(K)sABS(CV1(K)) 9 ABSOLUTE RELATIVE INLET VELOCm (6) 
AV1(R):PHASE(CV1(K))-i>90 9 ANGLE OF RELATIVE INLET VEL (7) 

§••§• impingement angles ••••• 


AI0=A1(K)4-AV1(K)-90 9 INLET IMPINGQfENT ANGLE (BELOW SURFACE) 

IF (AIO.GE.O) THEN 

IF (NOPT.EQ.O) PRINT IMP*,AI0,' BELOW SURFACE IN STG»,K 
IF (NOPT.EQ.O) GO TO 10 9 ASK FOR NEW INPUTS 
ETAR=0 9 LOW VALUE TO TELL SEARCH THIS IS A BAD ANGLE 
GO TO 83 9 RETURN TO SEARCH DRIVER 
END IF 

NSTEP1=NSTEP+1 

AXUAI(K) 9 INITIAL GUESS AT ANGLE FROM CNTR BACK TO IMP PT 

DO 55 J^I.NSTEPI 9 FOR EACH STREAM TUBE PLUS LAST STREAMLINE 
YST=SB»(J-0.5)/NSTEP 9 CENTER DISTANCE OF STREAM TUBE 
IF (J.EQ.NSTEP1) YSTsSB 9 LAST STREAMLINE 
C1=XSIN(A1(K))+(XC0S(A1(K))-YST/R1(K))»XTAN(AV1(K)) 9 (10) 

DO 50 JJ=1,100 9 ITERATE EQ.(9) UP TO 100 TIMES FOR AX 
AX=PXX*ASIN(C1 -XTAN(AV1(K))«XC0S(AX1)) 9 NEW AX (9) 

IF (ABS(AX-AXI).LT.O.OI) GO TO 55 9 CONV TO 0.01 DEG 
CALL RGLFSI (AX1,AX,JJ) 9 CALL REGULA-FALSI ITERATION 

PRINT •,» DID NOT CONVERGE ON AX IN STAGE’ ,K 9 AFTER 100 TIFfiS 

THETA(J)=90-AV1(K)-AX 9 IMPINGQIENT ANGLE (11) 

IMPACT AND FRICTION LOSS ••••• 

V0=0 9 DUMMY VALUE OF INLET VELOCITY FOR FIRST SECTION 
C 

DO 60 J=1,NSTEP 9 CALCULATE VELOCITY LEAVING EACH SECTION 
M0=(J-1)»DM 9 FLOW RATE AT INLET OF SECTION 
MM=M0+DM/2 9 FLOW RATE AT MIDDLE OF SECTION 
ME=J»DM 9 FLOW RATE AT END OF SECTION 
C 

REM=1*»MM/(HN»VISCL) 9 REYNOLDS NUMBER (17) 

CALL FRIC(REM.CFM) 9 CALCULATE FRICTION COEFFICIENT 
C 

Z1=( 1.0 + M0»V0/(DM»V1(K)»XC0S(THETA(J))) ) •(DM/ME)»»2 
Z2=XSIN(THETA(J))»ME/(CFM»RAP1»DM) 
TERM2=2»CFM»RAP1»Z1/XTAN(THETA(J)) 9 2ND TERM IN SQRT 
IF (TERM2.lt. IE-12) THEN 
IF (NOPT.EQ.O) PRINT TERM2 TOO SMALL’ 

IF (NOPT.EQ.O) GO TO 10 9 ASK FOR NEW INPUTS 
ETARsO 9 LOW VALUE TO TELL SEARCH THIS IS A BAD ANGLE 
GO TO 83 9 RETURN TO SEARCH DRIVER 
END IF 
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c 

c 

60 

C 

C 

C 

C 


c 


c 


c 

c 

c 


c 

c 

c 


c 


c 


VE=V1(K)»Z2»( DSQRT(1.0D0+TERM2) -1.0DO) t SECTION EXIT VEL (28) 

IF (NDET.GE.7) PRINT 00T6 

VOsVE 8 INLET VELOCITY FOR NEXT SECTION 

V2(K)rVE 8 LIQUID VELOCITY LEAVING IMPINGQ1ENT ZONE 

••••• FRICTION LOSS IN FILM FLOW ZONE ••••• 

REF(K):4*HB/(HN*VISCL) 8 REYNOLDS NUMBER IN FILM ZONE (30) 

CALL FRIC( REF(K),CFF(K) ) 8 FRICTION COEFF IN FILM ZONE 

IF(N0PT.EQ.0.AND.-AX.GT.A3(K)) PRINT*, » IMP BEYOND Rl, STG'.K 
IF (N0PT.GT.1.0R.NR1.EQ.1) A3(K)=-AX 8 SET END OF Rl AT LAST STRMLN 
L1(K)=R1(K)*(A1(K)+A3(K))*PX 8 LENGTH TO END OF INLET RADIOS (3D 
LIMP(K)=R1(K)*(A1(K)-AX)*PX 8 LENGTH TO END OF IMPINGEMENT ZONE (32) 
LS0RF(K)=L1(K)+R2(K)*(A2(K)-A3(K))*PX+LEXT(K) 8 TTL SURF LENGTH (33) 
LFILM=LSURF(K)-LIMP(K) 8 LENGTH OF FILM-FLOW ZONE (34) 

V3(K)=1/( 1/V2(K)+0.5*RH0L*CFF(K)*HN*LFILM/MB ) 8 REL EXIT VEL (36) 
TFILM=MB/(RH0L*HN*V3(K)) 8 FILM THICKNESS LEAVING BLADE 
V3(K)=V3(K)»XCOS(ADIV) 8 CORR FOR FLOW DIVERGENCE IF ADIV>0 (57) 

ttttfft blade EXIT CONDITIONS ••••• 

CV3(K)=POLAR(V3(K),A2(K)) 8 VCTR RELATIVE EXIT VELOCITY (37) 
CVOOT(K)=CV3(K)+CVB(K) 8 VCTR ABSOLUTE EXIT VELOCITY (38) 
VOOT(K)=ABS(CVOOT(K)) 8 ABSOLUTE EXIT VELOCITY (39) 
AOUT(K)=90+PHASE(CVOOT(K)) 8 ANGLE OF ABSOLUTE EXIT VELOCITY (40) 

••••• WINDAGE TORQUE ••••• 


IF (K.EQ.D WREL1:VB(K)/RB 8 REL INLET ANG VEL FIRST STAGE 
IF (K.GT.D WREL1=(VB(K)-VB(K-1))/RB 8 STAGES AFTER FIRST 
IF (K.LT.NSTG) WREL2=( VB(K)-VB(Kt-1 ) )/RB 8 STAGES BEFORE LAST 
IF (K.EQ.NSTG) WREL2 j^VB(K)/RB 8 REL EXIT ANG VEL LAST STAGE 

RE1=RH0G*ABS(WREL1)«RW«»2 /VISCG 8 REYNOLDS NO., INLET SIDE (41) 
RE2=RH0G*ABS(WREL2)«RW»»2 /VISCG 8 REYNOLDS NO., EXIT SIDE (41) 
RE1:MAX(RE1,100) 8 TO PREVENT OVERFLOW IN CM01 AT ZERO SPEED 
RE2=MAX(RE2,100) 8 TO PREVENT OVERFLOW IN CM02 AT ZERO SPEED 
CC=R1(K)*( XSIN(A1(K))+XSIN(A3(K)) ) ♦ R2(K)«( XSIN(A2(K)) 

& - XSIN(A3(K)) ) LEXT(K)*XC0S(A2(K)) 8 BLADE CHORD (42) 

AR=HB/CC 8 ASPECT RATIO OF BLADES 

CMSTARsO. 0067+0. 007*AR»»0. 6 8 FIT TO MANN & MARSTON TABLE 3 (43) 
GAMMA=MAX(0. 21-0. 086«AR, 0.15) 8FIT TO MANN & MARSTON TABLE 3 (44) 

CM01=CMSTAR»((1E6/RE1)**GAMMA)/(1+2.3*0.14) 8 AT RE1 , CCaO (46) 
CM02=CMSTAR»((1E6/RE2)*»GAMMA)/(1+2.3*0.14) 8 AT RE2, CC^^O (46) 
CM1sCM01*(1+2.3*CC/RW) 8 BLADED MOMENT COEFF, INLET SIDE (45) 
CM2=CM02»(1+2.3*CC/RH) 8 BLADED MOMENT COEFF, EXIT SIDE ^45) 
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c 

M1«SIGN(1.0,WREL1)»0.25»CM1*RH0G»WHEL1»»2 "RW^^S « INLET HINDG (47) 
M2«SIGN(1.0.WREL2)»0.25»CM2»RHOG»WREL2»»2 •RW^^S « EXIT WINDG (47) 

C 

LW(K)sMUM2 0 TOTAL WINDAGE TORQUE ON STAGE K 
C 

C ••••• LIQUID TORQUE ••••• 

C 

FL(K)sML»(VIN(K)»XCOS(AIN(K))-(1-FSTAG)»VOUT(lO»XCOS(AOUT(K)) 

A -FSTAG«VB(K)) 9 LIQUID FORCE (59) 

C 

LL(K)sFL(K)«RB 9 LIQUID TORQUE (49) 

70 IF (K.LT.NSTG) CVIN(K+1 )=(1-FSTAG)»CVOUT(K) 

& ♦ FSTAG*CVB(K) § VECTOR ABSOLUTE INLET VEL TO NEXT STAGE (60) 

C 

c ••••• gas-phase torque *••*• 

C 

DO 75 Ksl.NSTG 0 CALCULATE GAS- PHASE TORQUE 
IF (K.EQ.1) VGK=VG 0 INLET GAS VELOCITY TO FIRST STAGE 
IF (K.GT.1) VGKsVIN(K) 0 USE GAS VEL EQUAL LIQ VEL AFTER FIRST STAGE 
LG(K)s2»ETAG»MG»VGK»RB»(1.0-VB(K)/VGK) 0 GAS-PHASE TORQUE (50) 
LB(K)sLL(K)4LG(K) 0 NET BLADE TORQUE 
75 LR(K):LB(K)-LU(K) 0 NET ROTOR TORQUE 

C 

C §•§•• POWER ••••• 

C 

PB=0.0 0 INITIALIZE TOTAL BLADE POWER TO ZERO 
PRsO.O 0 INITIALIZE TOTAL ROTOR POWER TO ZERO 
C 

DO 80 K=1.NSTG 0 SUM POWERS OVER STAGES 
PBK=LB(K)*WW(K) 0 BLADE POWER IN STAGE K (51) 

PRK=LR(K)*WW(K) 0 ROTOR POWER IN STAGE K (52) 

PB=PB+PBK 0 TOTAL BLADE POWER 
80 PR=PR+PRK 0 TOTAL ROTOR POWER 

C 

C ••••• EFFICIENCY ••••• 

C 

ETAB=PB/PJET 0 BLADE EFFICIENCY (54) 

ETAR=PR/PJET 0 ROTOR EFFICIENCY (55) 

ETAT=ETAR»ETAN 0 TURBINE EFFICIENCY (56) 

C 

83 IF (NMON.EQ.1) PRINT 85.ETAR 

IF (NMON.EQ.2) PRINT 85,ETAR, (X( J) , J=1 ,NX) 

85 FORMAT (F1 1 .8,6(F10.5) ) 

Y=1-ETAR 0 QUANTITY TO BE MINIMIZED IN SEARCH 
IF (NOPT.GT.O) GO TO 30 0 FORM NEW X-VECTOR 
C 

90 IF (KGO.GT.2) PRINT KGO=',KGO 0 ERROR 
C 

PRINT OUT1 
C 

IF (NDET.GE.2) PRINT 0UT2 
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IF (NDET.GE.3) PRINT 00T3 
IF (NDET.GE.H) PRINT 0UT4 
IF (NDET.GE.5) PRINT 00T5 
C 

95 GO TO 10 9 ASK FOR NEW INPUTS 
C 

100 STOP 


•fftitititt FRICTION COEFFICIENT SUBROUTINE •#•§§••#§ 
SUBROUTINE FRIC(RE.CF) 

IF (RE. LT. 1034.8) CF:16/RE § LAMINAR FLOW (18) 

IF (RE.LT.1034.8) RETURN (f LAMINAR FLOW 

CFGS=0.001525+0.1375/CBRT(RE) « APPROXIMATE CF 

DO 10 KK=1,100 

CF=1/(4»L0G10(2»RE»SQRT(CFGS))-1.6)»»2 § VON KARMAN EQ. (19) 
IF (ABS(CF-CFGS).LT.IE-T) RETURN # CONVERGED 
0 CFGSrCF e NEXT CF GUESS 

PRINT • DID NOT CONVERGE ON CF* 

RETURN 


REGULA-FALSI ITERATION SUBROUTINE 


C 


C 


C 

C 

10 

15 

C 

C 


SUBROUTINE RGLFSKXI ,X2,K0UNT) 

IF (KOUNT.EQ.1) SL0PE=0.0 § SLOPE FOR FIRST PASS 

IF (KOUNT.GT.1) SL0PE=(X2-X20LD)/(X1-X10LD) § SLOPE FROM PREV POINT 

X10LD=X1 8 NEW ABSCISSA BECOMES OLD ABSCISSA 
X20LD=X2 # NEW ORDINATE BECOMES OLD ORDINATE 

X1=(X2-SL0PE«X1)/(1D0-SL0PE) € ABSCISSA CROSSING 45-DEG LINE 

IF (NM0N.EQ.3.AND.K0UNT.EQ.1) PRINT 10 
FORMAT(/5X, *X1 ' ,11X, 'X2' ,9X, 'SLOPE'/) 

IF (NMON.EQ.3) PRINT 15, XI ,X2, SLOPE 
F0RMAT(3(E13.6)) 

RETURN 

END 


125 


Table A-1. Rotor model nomenclature 


Program 

Text 

Definition 

A1 


blade inlet angle (Fig. A-1) 

A2 


blade exit angle (Fig. A-1) 

A3 


angle to end of inlet radius (Fig. A-1) 

ADIV 

div 

average divergence angle of exit flow 

AIO 


impingement angle at blade inlet 

AIN 

A. 

in 

angle of absolute inlet velocity (Fig. A-2) 

ANOZ 

A 

noz 

nozzle angle (Fig. A-1) 

AOUT 

A ^ 
out 

angle of absolute exit velocity (Fig. A-2) 

AR 

AR 

blade aspect ratio H^/C 

AVI 

Ayi 

angle of relative inlet velocity (Figs. A-1, A-2) 

AX 

A 

X 

angle to impingement point (Fig. A-3) 

AXl 


guess at AX 

Cl 

S 

constan." term in Eq. (A-9) 

CC 

c 

blade chord (Fig. A-3) 

CF 


friction coefficient 

CFF 


friction coefficient in film-flow zone 

CFGS 


guess at 

CFM 

^fm 

friction coefficient in middle of blade section 

CMOl 

Sn 

0 

windage moment coefficient at zero thickness, 
inlet side 

CM02 

C„, 

0 

windage moment coefficient at zero thickness, 
exit side 
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Program 


Text 


Definition 


CMl 

Cm 

windage moment coefficient on inlet aide 

CM2 

Cm 

windage moment coefficient on exit aide 

CMSTAR 

c ^ 

windage moment coefficient for * 10^ 

and C/R„ - 0.14 

CVl 

Vl 

vector relative inlet velocity (Fig, A-2) 

CV3 

V 3 

vector relative exit velocity (Fig, A-2) 

CVB 

Vb 

vector blade velocity (Fig, A-2) 

CVIN 

Vin 

vector absolute inlet velocity (Fig. A-2) 

CVOUT 

Vout 

vector abaolute exit velocity (Fig. A-2) 

DM 

dm 

liquid flow rate impinging on blade section 
(Fig. A-5) 


Dh 

hydraulic diameter 

EPSF 


convergence criterion on dependent variable 
in search routine 

EPSX 


convergence criterion on independent variable 
in search routine 

ETAB 

’’b 

blade efficiency Pb/^jet 

ETAG 

"g 

gas torque factor (peak gas-phase efficiency) 

ETAN 

'’n 

nozzle efficiency 

ETAR 

'’r 

rotor efficiency Pf/Pjet 

ETAT 

'’t 

turbine efficiency 

FL 

Ft 

force of liquid phase 

FSTAG 

^stag 

specified fraction of liquid flow stagnated 
in the rotor 
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Table A-1 (conCd) 


Program 

Text 

Definition 

GAMMA 

y 

exponent of Reynolds number in moment 
coefficient 

HB 

Hb 

blade height (Fig. A-1) 

HN 

Hn 

nozzle height (Fig. A-1) 

ITMAX 


maximum number of iterations in search routine 

KGO 


flag used to control search routine 

Ll 

h 

length to end of inlet radius (Fig. A-6) 

LB 

s, 

blade torque 

LEXT 

^ext 

length of straight extension (Fig. A-1) 

LFILM 

ilm 

length of film flow zone (Fig. A-6) 

LG 

L 

g 

gas-phase torque 

LIMP 

L. 

imp 

length of impingement zone (Fig. A-6) 

LL 


liquid-phase torque 

LR 

L 

r 

rotor torque 

LSURF 

^surf 

blade surface length (Fig. A-6) 

LW 

L 

w 

windage torque 

MO 


liquid flow rate entering blade section 
(Fig. A-5) 

Ml 


windage torque on inlet side 

M2 

”2 

windage torque on exit side 

MB 


liquid flow rate on each blade 

ME 

“e 

liquid flow rate leaving blade section 
(Fig. A-5) 
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Table A-1 (contd) 


Program 

Text 

Definition 

MG 

®g 

gas flow rate 

ML 

“t 

liquid flow rate 

MM 


liquid flow rate in middle of blade aection 
(Fig. A-5) 

MXSTEP 


maximum number of steps for blade impact and 
friction calculation 

NDET 


selects level of detail of printout 

NINPT 


1 for printout of inputs 

NMON 


1 or 2 for monitoring of search 

NOPT 


1 for rpm optimization, 2 for rpm and blade 
angle optimization 

NRl 


selects matching of blade inlet length to 
impingement length 

NSTEP 


number of blaae steps for impact and friction 
calculat ion 

NSTEPl 


NSTEP ♦ 1 

NSTG 


number of rotor stages 

NX 


number of variables to be optimized 

PB 


total blade power 

PBK 


blade power in stage K 

PI 

ir 

pi 

PJET 

^ jet 

2 . 2 

jet power S (mV. ♦mV ) 

PR 


total rotor power L^u> 

PRK 

Pr 

rotor power in stage K 




Table A*1 (contd) 


Program 

Text 

Definition 

PX 


*/180 

PXX 


180/» 

R1 

“l 

blade inlet radius (Fig. A-1 ) 

R2 


blade exit radius (Fig. A*l) 


R 

a 

gas/liquid flow area ratio 

RAPl 


R ♦ 1 
a 

RB 

Rb 

rotor radius at center of nozzle exit 
(Fig. A-l) 

RE 

R* 

Reynolds number 

REl 

Rej 

Reynolds number of windage flow on inlet side 

r.E2 

R*2 

Reynolds number of windage flow on exit side 

REF 

R*f 

Reynolds number of film flow 

KEM 

Re 

fn 

Reynolds number in middle of blade section 

RHOC 

0 

g 

gas density 

RHOL 

‘’t 

liquid density 

RPM 


rotor speed, revolutions per minute 

RW 

R 

w 

outside radius of rotor (Fig. A-l) 

SR 


blade spacing (Fig. A-l; 

TERM2 


second term under the square root in 
Eq. (A-28) 

TFILM 


thickness of liquid film leaving blade 


^m 

thickness of liquid film in middle of blade 
section (Fig. A-5) 

THETA 

6 

impingement angle (Fig. A-3) 
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Table A-1 (contd) 


Program 

Text 

Definition 

VO 

Vo 

velocity of liquid film entering blade section 
(Fig. A-5) 

VI 

Vl 

relative inlet velocity (Fig. A-1) 

V2 

V2 

velocity of liquid leaving impingement zone 
(Fig. A-1) 

V3 

V, 

relative exit velocity (Fig. A-1) 

VB 

Vb 

blade velocity (Fig. A-1) 

VE 

Ve 

veloc i: liqaid leaving blade section 

(Fig. A 

VG 

Vg 

gas veloc ic 

VGK 


inlet gas velocity to stage K 

VIN 

^in 

absolute inlet velocity (Fig. A-2) 


Vm 

velocity of liquid in middle of blade section 
(Fig. A-5) 

VISCG 


gas viscosity 

VISCL 


liquid viscosity 

VL 

Vl 

liquid velocity 

VOUT 

Vout 

absolute exit velocity (Fig. A-2) 

WA 


array used in search routine 

WN 

Wn 

nozzle width (Fig. A-1) 

WRELl 


speed of rotor relative to previous stage or 
atmosphere 

WREL2 


speed of rotor relative to next stage or 
atmosphere 
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Table A-1 (contd) 


Program Text 


Definition 


ww 

(D 

rotor speed, radians per second 

X 


array used in search routine 

XBEST 


array used in search routine 

Y 


1 - n^, quantity minimized in optimization 
search 

YBEST 


variable used in search 

YST 

Yst 

distance of streamline from blade inlet 
(Fig. A-3) 

Z1 


product of last two factors in second term 
under the square root in Eq. (A-28) 


Z2 


product of first two factors in Eq. (A-28) 





APPENDIX B. 


SEPARATOR TURBINE MODEL 


Figure B-1 shows the geometry of the separator turbine. The two-phase 

flow enters at radius R and angle A • The liquid has flow rate 

noz ° noz 

m and velocity The gas has flow rate m and velocity V • 

The jet power is 


P. 

jet 


ifm V ^ + m V 
\ I t g g / 


(B-1 ) 


A. Separator Flow 

The liquid impinges on the inside wall of a rotary separator turning at 
angular velocity The liquid forms a spinning liquid layer of radius 

Rsep’ liquid is removed from the separator by a liquid turbine with one 

or more U-tube scoops turning at angular velocity The speed of the 

liquid layer is 


V = uj,R (B-2) 

sep 1 sep 


B. Liquid Turbine Flow 

The speed of the scoop is 


V ^ = o)^R 
turb 2 sep 

The relative velocity of the liquid entering the scoop is 

V = V - V 
rel sep turb 


(B-3) 


(B-4) 


The relative velocity of the liquid leaving the scoop is • The 

scoop efficiency can be defined as the square of the ratio of exi t ^ve 1 oc i ty to 
inlet velocity. Thus, if the scoop efficiency is given, the 

relative exit velocity is 
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Fig. B-1. Separator turbine 


134 



(B-5) 


V , * V , x/n 

rel^ rel^\ scoop 

The absolute forward exit velocity of the liquid leaving the liquid 
turbine is 


V * V - V 
out turb rel 2 


The required inlet area of the scoop, or scoops, is 


(B-6) 


A 


scoop 


ni. 


V , 0, 
rel^ 1 


(B-7) 


where is the liquid density. 

The torque exerted on the scoop by external liquid drag can be expressed 
as a drag coefficient times the product of frontal area and dynamic 
pressure : 


ext 


bC.A p V , R 
d scoop I rel^ sep 


(B-8) 


C. Liquid Torque 


The angular momentum of the liquid leaving the nozzle is 


M, = m V R cos A 
1 I i noz noz 


(B-9) 


The angular momentum of the liquid leaving the separator is 


* mV R 
2 1 sep sep 


(B-10) 


The angular momentum of the liquid leaving the liquid turbine is 


M » m V R 
3 I out sep 


IIS 


(B-11) 




The sepatator torque due to the liquid flow is equal to the change in 
angular momentuin of the liquid in the separator minus the external drag of the 
liquid turbine: 


h 


iq 


M - - L , 

1 2 ext 


(B-12) 


The liquid turbine torque is equal to the change in angular momentum of 
the liquid in the turbine, plus the external drag. 


L , = L 

turb 2 3 ext 


(B-13) 


D. Windage 


The windage torque of the separator can be calculated from the windage 
torque formulas of Appendix A, evaluated for a nonbladed disc. The Reynolds 
number is 


Re 


%“l"o 


(B-14) 


where 

u is 
g 


0 ^ is the gas density, 
the gas viscosity. 


R is 
o 


the separator outside radius. 


and 


With no blades, the aspect ratio AR is zero. From Eq . (A-A3), the moment 
coefficient for a disc of thickness/radius ratio = 0.14 at Re = 10^ is 
= 0.0067. From Eq . (A-44), the exponent of the Reynolds number 
dependence is y = 0.21. From Eq . (A-46), the moment coefficient for a thin 
disc at Re = 10^ is then 


C 


m 


0 


0.0067 (lO^/Re)*^'^^ 
1 + (2.3H0.14) 


(B-15) 
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From Eq. (A-45), the moment coefficient for the separator of width 
is 

C “ C f 1 + 2.3 I (B-16) 

•"o "b \ ^0 / 

From Eq. (A-47), the windage torque is 

L = 'i C 0 11?, (B-17) 

w m g I 0 

The windage torque of the liquid turbine is ignored because of its lower 
speed • 


E. Gas Torque 


The efficiency of power recovery from the gas phase at optimum gas 
turbine speed is a specified constant, If a separate gas turbine is 

used that runs at optimum speed, then the gas turbine power is 


P 


gt 


^ n m V 
g g g 


(B-18) 


If the gas turbine consists of impulse blading attached to the separator 

at radius R , then n is the fraction of ideal gas torque at the 
gt g 

separator speed. The gas torque is thus 


L 

gas 


2 n m V R ^ cos A 
g g g gt noz 



0) 

sync 


U-19) 


where ^ is the synchronous speed at which the blade speed equals the 
circumferential component of gas velocity. 


U) 

sync 


V cos A 
g no z 


(B-20) 
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F 


Power and Efficiency 


The total separator torque is 


L - L . ♦ L - L 

sep iiq gas w 


(B-21) 


where L is zero if a separate gas turbine is used, 
gas 

In the calculation presented in Fig, 43, L was zero and the separator 

ga 8 

speed was adjusted to make L zero (free-spinning separator). 

sep 


The separator power output is 


P “ L (D, 
sep sep 1 

The liquid turbine power output is 

^turb ^turb ^2 

The total power output is 


(B-22) 


(B-23) 


where P ^ 

gt 

Th e 


P = P P . + P 

tot sep turb gt 


is zero if the gas turbine is attached to the separator. 


rotor efficiency is 


(B-24) 


and the turbine efficiency is 


tot 

^jet 




where is the nozzle efficiency. 


(B-25) 


(B-26) 
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